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INTRODUCTION 

This  Workshop  is  intended  as  a  small  meeting  of  people  active  in  the  field 
of  the  interaction  of  hydrogen  with  dopants  and  defects  in  III-V  compounds  and 
their  alloys.  New  results  and  analyses  will  be  presented  and  discussed  and  it  is 
hoped  this  workshop  will  provide  a  better  understanding  of  the  different  ways 
hydrogen  can  be  introduced  in  III-V  semiconductors,  of  the  phenomena  related 
to  its  presence  and  of  its  possible  use  to  improve  the  quality  of  devices. 

The  main  subjects  debated  in  this  workshop  will  be: 

-  The  practical  use  of  hydrogenation,  either  by  plasma  introduction  or  bv 
proton  implantation  in  the  fabrication  of  III-V  devices. 

-  The  study  by  electrical  and  optical  methods  of  the  neutralization  bv 
hydrogen  of  dopants  in  III-V  compounds. 

-  The  correlation  between  the  distribution  of  hydrogen,  its  charge  state  and 
the  dopant  neutralization. 

-  The  nature  and  properties  of  point  defects  and  extended  defects  related  to 
proton  implantation. 

-  The  consequences  of  the  presence  of  hydrogen  in  bulk  ,  MOVPE  and  MBE 
materials  non  intentionally  hydrogenated. 

The  presentation  of  channeling  measurements  which  give  a  good  insight  of 
the  atomic  location  of  hydrogen  in  the  crystal  lattice  is  a  natural  transition  to  the 
calculation  of  microscopic  models  of  the  hydrogen  complexes,  first  performed 
for  hydrogen  in  pure  and  doped  silicon  and  more  recently  in  doped  gallium 
arsenide. 

We  wanted  finally  this  workshop  to  include  a  presentation  of  the  new 
important  results  on  muonium  level  crossing  spectroscopy  in  silicon,  gallium 
arsenide  and  gallium  phosphide.  They  provide  information  on  the 
configuration  of  anomalous  muonium  in  these  materials  and  the  link  with  the 
recent  EPR  results  on  isolated  hydrogen  in  silicon  could  be  considered  as  a 
serious  experimental  indication  of  the  bond-centred  nature  of  hydrogen  in  the 
pure  materials. 

This  booklet  contains  the  abstracts  or  extended  abstracts  of  the  talks  by  the 
invited  speakers  and  the  reports  we  received  from  the  other  participants.  We 
have  tried  to  classify  them  as  follows: 

1)  Applications  and  technology. 

2)  General  studies  and  electrical  measurements. 

3)  Extended  and  point  defects. 

4)  Vibrational  spectroscopy. 

5)  Advanced  nuclear  methods. 

6)  Theory. 


HYDROGEN  PASSIVATION  OF  DOPANTS  AND 
DEFECTS  IN  lll-V  COMPOUNDS  AND  THEIR  ALLOYS 


Applications  and  Technology 


N.  Caglio,  E.  Constant  (Universite  de  Lille),  J.  Chevallier  and  J.  C.  Pesant:  (C.N.R.S. 
Bellevue)  "A  new  process  for  the  fabrication  of  field  effect  transistors  using 
the  neutralization  of  shallow  donors  by  atomic  hydrogen". 

N.  Duhamel  (C.N.E.T.  Bagneux):  "Electrical  insulation  by  proton  bombardment 
in  III-V  materials:  Application  to  micro-  and  opto-electronics". 

S.  Cole,  L.  \Y.  J.  Duncaii,  W.  j.  M.  Rcthuvii,  P.  J.  Skcvir.gtor.  a;.d  G.  T.  D. 

Spiller  (British  Telecom  Ipswich):  "The  effect  of  cooling  ambient  on  the 
electrical  activation  of  dopants  in  MOVPE  of  InP". 

G.  R.  Ante!!,  A.  T.  R.  Briggs,  B.  R.  Butler,  S.  A.  Kitching,  J.  P.  Stagg  (STC 
Technology  Harlow),  A.  Chew  and  D.  E.  Sykes  (University  of  Tediuology 
Loughborough):  "The  passivation  of  Zn  acceptors  by  hydrogen  during 
MOVPE  growth". 

Ph.  Collot  and  Ch.  Gaonach  (Thomson  CSF  -  LCR  Orsay):  "Non-intentional 
neutralization  of  dopant  in  GaAs  by  CH4/H2  reactive  ion  etching". 
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:nthcduct::n 


'■.hen  atomic  hydrogen  is  diffused  into  .-.-type  silicon  doped  Z 
fcrr.at ion  of  silicon  related  hydrogen  complexes  occurs  [.'alii  S.  al .  13 
As  a  result,  the  concentrator  of  active  silicon  donors  drops  do-.-,  as 
through  the  decrease  of  the  free  carler  concer.trat ion  '-hich  is  acccm.ca 
by  a  simultaneous  increase  of  the  electron  nobility.  It  is  the  goal  of 
paper  to  describe  how  to  realize  field  effect  transistor  by  using 
phenomena.  In  a  first  part,  the  fabrication  process  first  used  to  rea 
CaAs  MHSFF7  is  presented,  then,  we  describe  the  method  of  modeling  w 
has  been  used  to  optimize  this  technological  process,  lastly,  we  she 
descrioe  tne  first  realization  of  AlCaAs/CaAs  MI  Sr  £  I  based  upon 
neutralisation  cf  donors  by  atomic  hydrogen, 

b.  SFoT 


The  process  can  be  described  as  fellows,  '-'e  start  from  a  si 
highly  doped  CaAs  epi  layer.  After  the  mesa  ar.d  the  source  and  drain  c 
contacts  fabrication,  a  protecting  dielectric  layer  is  deposited  an: 
proceed  to  the  gate  lithography.  The  processed  sample  is  then  exposed 
hydrogen  plasma  for  hydrogen  diffusion  in  the  ncr.-protected  regions, 
hydrogen  diffusion  results  in  the  neutralisation  cf  a  mere  or  less  1 
fraction  cf  the  silicon  deners  depending  cn  the  plasma  exposure  paramei 
As  the  neutralisation  becomes  mere  effective,  the  electron  met; 
Increases.  After  the  gate  deposition,  the  devices  are  expected  to  have 


-4- 


Their  character! st :  cs  are  presented  cr.  Fig. 1.  The  gate  .er.gth  -as 
determi red  bv  ‘he  coer.irg  width  made  in  the  pciyimioe  protentton  -ayer.  it 


-■as  i .  2  pm  ah'  ter  exposure.  The  processed  sample  -as  tr.er.  -eft  :  r.  a  ?.,  ~ 
capacitive!/  coupled  hydrogen  plasma  for  -.3  r.ir.utes  at  ZIZ'Z  ;  in  such 
conditions,  only  a  few  per  cent  cf  the  total  co.nror  concentration  re-air.s 
active  after  exposure,  -e  then  proceeded  to  a  snail  recess  of  3.  1  pm  ir.  the 
gate  region  before  the  gate  deposition.  Three  different  geometries  have 
been  used  with  gate  width  cf  75,  150  ar.d  300  pm.  We  finally  annealed  the 
transistor  at  300'C  for  15  minutes  in  order  to  adjust  the  active  do.rtr 
concentration  in  the  active  region  to  the  desired  value.  .Actually,  it  has 
been  shown  that  the  .neutralisation  effect  is  stable  for  several  years  or 
mcre  at  130'C  ar.d  below,  ar.d  that,  above  2S0*C  annealing  leads  to  a 

deneut  rai  i  sat  ion  of  the  doncns.  A  very  simular  process  has  been  -sed  to 
reiize  0.4  pm  gate  KFET  but  electronic  lithography  has  been  employed. 

figure  2  shews  the  static  I-V  characteristics  of  the  HrET  with,  a  310  ,-n 
gate.  The  corresponding  gm  value  measured  for  a  constant  so-p-ce-crai  n 
voltage  of  2  volts  is  as  high  as  2S0  xS/nm.  A  more  co-piete 

characterization  cf  these  devices  has  been  performed  by  nagme- ores : stance 

experiments  [Sites  &  al.  iSSOl.  Applying  this  technic  to  our  K FET ,  -e  find 
an  electron  mobility  cf  3500  cm  /v.s.  in  the  active  layer.  Taxing  into 
account  a  doping  level  cf  5.  13'  cm  3  deduced  f-cm  gate-source  capac i tan.ee . 
the  electron  mobility  value  is  satisfactorily  good.  The  equivalent  oirou.t 
cf  these  transistors  in  the  high  frequency  regime  has  been  deduced 
[Danbri.ne  !.  al,  1337).  The  cut-cff  frequency  is  higher  than  15  GHz  ar.d  35 
GHz  for  a  gate  length  of  1.2  urn  and  0.4  pm  respectively.  The  high  frequency 
tnansccnductar.ee  in  both  cases  is  higher  than  SCO  mS/m-  The  frequency 
dependences  of  the  Current  Gain.  Maximum  Au-qilzbie  Gain,  a  no  un:  lateral 
Gain  are  gi.en  cr.  figure  3.  for  a  i.Tp.r  gate  length  transistor,  the 

perfor-ances  appear  somewhat  remarkable 
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As  neutral isa* ion  In  Si  doped  AlGaAs  is  effective,  we  realize  a 
r.etal-ir.sulator- sesiccnductor  field  effect  trar.sistor  (MISFED.  This  FE7  is 
usually  fabricated  cn  an  ur.dcped  (AlCaAs/CaAs)  heterostructure.  --here  a 
twc-dir.ensicnnal  gas  (2DEG)  is  induced  in  a  potential  well  at  the  interface 
by  the  application  cf  the  positive  gate  voltage.  We  start  with  a  highly 
doped  AlGaAs  epilayer  (for  good  ohmic  contacts  and  low  access  resistances') 
as  describe  Fig. 5.  Before  gate  deposition,  we  expose,  with  SI  N  as 
dielectric  layer,  the  active  region  to  an  hydrogen  plasma,  in  order  to 
strongly  decrease  the  donor  density  (with  a  rate>19C).  In  this  first 
realization,  we  succeed  to  observe  a  negative  resistance  (rig. S3,  which  is 
an  important  phenomenon  of  this  component.  But,  we  must  say.  that  static 
and  hyperfrequency  performances  are  disapcint  ir.g.  We  hope,  in  a  next 
realization,  improve  these  results  with  making  a  change  in  the  layer 
structure. 

CCNCUJSICN 


We  have  presented  new  types  cf  field  effect  transistor  based  cn  the 
neutral isat ion  of  shallow  donors  by  atomic  hydrogen  in  the  active  region, 
the  static  and  the  high  frequency  characteristics  can  be  ccr.sicered  as  very 
encouraging.  A  mere  extended  study  is  necessary  to  optimize  the  process  and 
to  evaluate  its  potential  use  in  microelectronics . 

details  cn  the  different  points  cf  this  paper, 
work  ( Ml Sm  and  modeling  methcdl,  in  the  oral 
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uis'ru2..i  :  Sher.atic  device  structure  cf  the  KFE7 

Figure  2  :  DC  voltage  character  1st i cs  cf  200  p.T.  gate  KFF7.  The  gate  voltage 
varies  from  (*  200  c.V)  to  (-  600  r.V) 

Tiffure_3  :  Variation  of  the  Current  gain.  Maximum  Available  Oair.  arc 
Unilateral  gain  versus  the  frequency  for  1.2  uni  gate  length  I-IFFT 
FiEure_4  :  Cor.parison  between  static  characteristics  obtained  with  the 
ncdelisation  and  experiments 

rlgure^5  :  Sher.atic  cross-section  of  MISFET  structure 
Figure_6  :  DC  voltage  characteristics  of  100  4m  gate  H-MISFET 
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The  successful  design  ar.d  fabrication  ci  III-V  integrated  circuits  requires 
a  high  degree  of  electrical  insulation  between  closely  spaced  active 
devices.  This  car.  be  dene  by  : 

.  selective  implantation  of  the  desired  dopant  Iona  through  a  suitable  mask. 

.  mesa  etching  of  the  active  layer 
.  ion  implantation  induced  damage 

However,  the  two  first  techniques  produce  an  insufficient  insulation  which, 
for  example  in  GaAs  MESFET  technology!  leads  to  variation  in  insulation 
voltage,  or  to  the  well  known  detrimental  backsetting  effect.  This 
relatively  low  degree  of  electrical  insulation  can  result  in  significant 
circuit  design  cor.traints  including  an  increase  in  the  minimum  spacing 
rules . 

In  the  third  case,  which  we  will  detail,  non-dopant  ions  are  implanted  into 
the  active  layer,  the  crystal  lattice  is  damaged,  trap  centers  are  created 
which  compensate  the  shallow  dopants.  The  quality  of  insulation  depend  cn 
several  parameters  :  the  material  cf  the  active  layer  (GaAs,  In?, 
Ir.GaAs...),  its  doping  level,  the  implanted  species  (H,  D,  B,  0)  ar.d  the 
implantation  and  annealing  conditions. 

Two  test  structures  can  be  used  to  measure  the  resistance  of  implantation 
induced  damage  layers  :  a  vertical  geometry  test  structure  employed  by 
several  workers  /l, 2/  and  a  lateral  geometry  which  reflects  the  structures 
of  technological  interest  /3,£/.  In  the  first  geometry,  since  the  bulk  cf 
doped  III-V  material  has  a  much  lower  resistance  chan  the  damage  layer, 
values  of  resistivity  can  easily  be  deduced  from  the  measure  of  the  current 
flow  frorr  the  front  to  the  back  contact  of  the  sample  as  a  function  c: 
applied  voltage.  In  the  second  configuration  only  insulation  re9istar.ee  car. 
be  deduced  without  any  hypothesis  cn  layer  thickness  and  current  path  but 
this  resistance  is  the  most  important  parameter  to  character ice  an 

insulation  for  mi ;r :elec trcr.ic  ccmponancs  . 

GaAs  In  GaAs,  the  insulation  resistance  is  higher  than  :: 

after  prot'.n  bombardment  /  5/  cr  Borer,  implantation  Thus  protons  or 

Hor?n  io“S  nave  been  successfully  employed  fst  device  insulation  in  GaAs 


.'CSFET  zizzziz*  /:,  3/.  However,  protons  have  the  disc  ^vantage  of  producing 
an  insulation  stable  ir.  a  snail  temperature  range  cr.ly.  Boron  icr.s  shew 
temperature  stability  up  to  500*C  / 1  * 6 / ,  while  protons  reach  only  -.OC^C, 
except  for  multi-energies  bombardments  /?/.  Another  point  to  be  underlined 
is  the  considerable  decrease  in  insulation  for  doses  higher  than  an  optimum, 
dose  due  to  an  enhanced  hopping  conduction.  In  the  case  of  active  layer 

17-  -  ;  1 4  «.*' 

MISFET  (r.»lxlO  e  cm  ),  the  optimum  dose  is  around  5x10  cm  "  for  protons 
i  1  -2 

/ 3 /  and  about  -xlO  -  n  for  Boron  ions  /£/.  This  means  that,  when  using 
Boron ,  the  implantation  time  is  reduced  by  three  orders  of  magnitude  for  the 
same  implantation  cur; ant. 

However,  cue  to  their  advantage  of  being  the  lightest  ions,  protons  are  used 
where  thick  damaged  layers  are  required.  For  example,  at  CNET  Laboratories, 
proton  bombardment  i«  routinely  performed  In  the  hetero June tion  bipolar 
transistors  (H5T)  technology  for  device  insulation  and  also  for  delimiting 
the  resistive  paths  into  the  n+  collector  layer  78/ .  It  has  also  been 
utilized  in  H3Ts  profess  to  reduce  extrinsic  base-collector  capacitance  and 
thus  tc  improve  their  '.igh  speed  performance  / 9 ,  10/. 

Finally,  ve  discuss  --ne  of  the  first  use  of  proton  bombardment  in  devices 
i.e.  the  formation  -.e  stripe-geometry  GaAs/GaAlAs  heterostructure  lasers 
711/  In  a  simple  technological  process. 

Ir.?  and  a:.pyg  ;  Coming  now  to  other  III-V  compounds,  especially 
l.nP  and  Ir.GaAs.  :cr  *he  purpose  of  interdevice  insulation,  implantation 
induced  damage  is  cv-nparable  in  GaAs  and  p  type  InP  ;  resistivities  of 
10-*. cm  can  be  ob  taint-!  by  proton  bombardment  / 1 2 / -  However,  in  n-type  Ir.?, 
the  maximum  resistivity  which  has  been  reported  for  proton  bombardment  is 
about  2x10  G.cn  /ll.Vnilt  in  In  -GalnAs  the  increase  of  sheet  resistance 
after  hydrogen  or  Bu*on  implantations,  is  too  small  to  produce  effective 
insulation  in  most  otical  or  electronic  InGaAs  devices  7127.  Thus,  the 
insulation  is  one  c*  the  major  problems  to  be  solved  before  a  viable 
technology  of  medium  complexity  circuits  will  be  established  ir.  these 
materials.  Thereby,  i  InP  technology,  protcn  bombardment  has  only  been  used 
ir.  phocrd* teccicn  71-  or  in  lasers  / 15/,  i.e.  in  p-cype  In?  material. 
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effect  on  surface  morphology. 


Samples  were  characterised  using  the  electrochemical  profiler  '  1  ’ 
ar.d  secondary  icn  mass  spectrometry  (SIMS)  .  Both  MCVPE  grcw.t  layers 
ar.d  Crochralski  bulk  crystals  were  studied  and  found  to  behave 
similarly.  It  was  found  that  if  AsH3  was  present  in  the  ceding 

ambient,  the  coping  level  was  lower  by  up  to  an  order  cf  magnitude 
than  if  ph3  was  used,  which  in  turn  cave  a  doping  level  lower  by  up 
to  a  factor  cf  2  than  obtained  with  K2  alone. r.-tvpe  In?  and  r.-  ;r 
p-type  GaAs  did  not  appear  to  be  affected. 

SIMS  analysis  showed  that  this  effect  was  caused  by  electrical 
deactivation  rather  than  loss  cf  the  dopant  species.  Significant 
quantities  cf  atomic  H  were  found  in  samples  exposed  to  the 
hydrides  during  cooling  (figure  1) ,  consistent  with  their  pyrolysis 
being  catalysed  by  the  semiconductor  surface.  Ke  suspect  these 
phenomena  to  be  caused  by  deactivation  of  the  dopant  by  complex 
formation  with  H,  as  observed  by  Johnson  et  ai.  [2]  in  SaAs  exposed 
to  a  hydrogen  plasma. 

On  introducing  p-dopant  precursors  into  the  cooling  ambient  a 
similar  lowering  of  doping  level  was  observed.  However,  SIMS 
analysis  could  detect  no  atomic  H  in  such  a  sample.  A  summary  cf 
some  of  our  findings  is  listed  below: 

1.  Doping  level  of  p-Inp  depends  on  cooling  ambient  such  that: 

p(AsH3)<p(PH3)<p(H2) 

2.  Doping  level  can  be  restored  by  heat-treating  in  ?H3  cr  H-2 

3.  Effects  common  to  MOVPE  epilayers  and  bulk  crystals 

4.  The  drep  is  not  due  to  loss  of  dopant 

5.  Magnitude  of  drop  is  proportional  to  original  doping  level 

6.  No  effect  seen  in  n-InP  or  n-  or  p-GaAs 

7.  Atomic  H  is  found  in  samples  cooled  in  hydride  ambients,  such 
that:  H(AsH3)>H(?H3) 

8.  Similar  effects  seen  in  samples  cooled  in  depar.t  precursor 
ambients 

9.  In  these  samples  the  original  dopant  diffuses  out  ar.d  is 
replaced  by  dopant  frem  the  arr.bie.nt-NO  H  IS  rC’JND 

To  conclude,  we  have  seen  that  the  electrical  activity  cf  p-dcpar.cs 
in  In?  is  strongly  influenced  by  the  cooling  ambient  and  that  the 
relative  concentrations  of  atomic  H  found  in  crystals  ceded  in 
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The  Passivation  of  Zn  Acceptors  in  In? 
by  Hydrogen  During  MOVPE  Growth 

G  K  Ar.teil,  ATS  Briggs.  B  R  Butler.  S  A  Kitching,  j  P  stagg. 
A  Chew*.  D  E  Sykes*  STC  Technology.  Harlow,  Essex. 
•Loughborough  Consultants  Ltd..  University  of  Technology 
Loughborough.  Leicestershire 


The  hole  concentrations  in  zinc  doped  InP  single  layers  grown 
by  MOVPE  are  reasonably  stable  when  annealed  at  about  400CC. 

When  a  similar  type  of  layer  is  grown  and  then  capped  by  a  layer 
of  lattice  matched  p  type  InGaAs  during  the  same  run  it  is  found 
that  the  hole  concentration  has  been  reduced  by  about  a  factor 
of  3  or  more  and  that  this  can  be  restored  to  the  higher  value 
by  annealing. 

At  first  it  was  thought  that  this  phenomenon  might  be  due  to  the 
rapid  diffusion  of  intenstitial  zinc  between  the  InP  and  InGaAs 
layers.  SIMS  analysis  of  a  specially  grown  structure  showed 
unequivocally  that  there  was  no  movement  of  zinc. 

Atomic  hydrogen  was  the  next  most  obvious  candidate  as  a  cause 
for  the  reduction  in  aCflfiBf.ftr  rv.rii’r-nl  i  *1  Inn  hiu  nan 
conci.io ted  to  t,«  vi.llKvJy  fur  the  roriowiny  reasons: 

1)  There  was  no  obvious  source  of  atomic  hydrogen  in  fhe 
MOVPE  growth. 

2)  The  InP  was  grown  in  a  predominantly  hydrogen  ambient  and 
yet  Inp  layers  could  be  grown  with  carrier  concentrations 
less  than  10'4cnr’  and  mobilities  in  excess  of 

200.000  volt  sec  cur’  at  77®K  and  thus  were  almost 
uncompensated.  To  produce  the  reduction  of  the  acceptor 
concentration  in  p  type  InP  would  require  mote  than 
lQl,cra~*  of  hydrogen. 

Despite  these  arguments  the  above  test  sample  was  re-examined  by 
SIMS  and  hydrogen  was  detected  in  the  p  type  InP  layers. 

Further  studies  soon  established  that  the  arsine  which  wa6 
present  during  the  growth  of  the  InGaAs  capping  layer  was  the 
source  of  the  atomic  hydrogen  and  that  this  hydrogen  was  trapped 
in  the  inP  mainly  during  the  cooling  down  stage  after  the 
completion  of  growth. 

It  Is  postulated  that  the  H-P  bond  in  InP  is  stronger  than  that 
of  the  As-H  bond  in  Astj  when  in  contact  with  the 
semiconductor  surface  and  thus  hydrogen  can  be  trapped  in  inP  at 
temperatures  where  Asf^  is  still  cracking. 

The  passivation  of  zinc  acceptors  in  InP  which  has  been 
overgrown  by  a  layer  of  p  type  InGaAs  can  be  avoided  by  the 
following  methods  : 
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1)  After  growth  is  complete  the  sample  can  be  cooled  as 
normal  in  an  AsH  ~Ft  ambient  to  about  500°C  at 
which  point  the  As^  is  switched  off  and  cooling 
continues  in  . 

2)  Grow  an  additional  layer  of  I nP  on  top  of  the  inGaAs  and 
cool  in  the  normal  Phj-Wj  ambient .  The 

additional  InP  layer  can  be  removed  by  a  selective  etch. 

3)  Grow  an  addtional  layer  of  n  type  InGaAs  on  top  of  the  p 
type  InGaAs  and  cool  as  normal  in  an  As/l 

ambient,  in  this  case  the  diffusion  rate  of  hydrogen  in 
n  type  InGaAs  is  much  slower  than  in  p  type  InGaAs  and 
thus  prevents  the  hydrogen  at  the  surface  from  teaching 
the  underlying  InP  layer. 
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HYDROGEN  PASSIVATION  OF  DOPANTS  AND 
DEFECTS  IN  lll-V  COMPOUNDS  AND  THEIR  ALLOYS 


General  studies  and  electrical  measurements 


W.  C.  Dautremont-Smith,  S.  J.  Pearton  and  M.  Stavola  (A.T.&tT.  Bell  Labs  Murray 
Hill):  "Recent  activities  in  hydrogenation  of  semiconductors".  Invited  talk. 

J.  Chevallier,  B.  Theys,  A.  Jalil,  R.  Rahbi  and  J.  C.  Pesant  (C.N.R.S.  Bellevue): 

Hydrogen  diffusion  and  dopant  neutralization  in  III-V  compounds 

Invited  talk. 

N.  M.  Johnson  (Xerox  Palo  Alto  and  Universitat  Erlangen-Numberg):  "Charge  states 
of  migrating  hydrogen  in  semiconductors".  Invited  talk. 

N.  Pan,  M.  S.  Feng,  K.  C.  Hsieh,  S.  S.  Bose,  G.  S.  Jackson,  G.  E.  Stillman  and  N. 
Holonyak,  Jr.  (University  of  Illinois  at  Urbana):  "  Impurity  and  structural  studies 
of  hydrogenated  GaAs  and  GaAs/Si". 

G.  Pensl  (Universitat  Erlangen-Numberg):  "Activities  on  hydrogen  passivation  of 
defects  in  silicon". 


RECENT  ACTIVITIES  LN  HYDROGENATION  OF  SEMICONDUCTORS 


\V.  C.  Dautremont-Smith,  S.  J.  Pearton  and  M.  Stavola 
AT&T  Bell  Laboratories,  Murray  Hill,  NJ  07974,  U.S.A. 


m-V  SEMICONDUCTORS 
Ga-As 

Passivation  of  donors: 

—  magnitude  and  depth  dependence  of  passivation,  dependence  on  donor  concentration 

—  thermal  stability/reactivation  kinetics,  dependence  on  donor  species,  proposal  of  donor  -  H 
direct  bonding  passivation  model 

—  H/D  in-diffusion  SIMS  depth  profiles,  dependence  on  donor  concentration,  effects  of 
various  doped  epi-overlayers 

—  variation  in  correlation  between  elemental  and  electrically  active  depth  profiles  with  type  of 
plasma  exposure  (low  frequency  vs.  high  frequency  vs.  microwave) 

—  defect  generation 

—  depth  dependence  of  electron  mobility  changes 

—  depth  dependence  of  donor  reactivation 

Passivation  of  acceptors: 

—  magnitude  and  depth  dependence  of  passivation,  efficiency  relative  to  donor  passivation 

—  thermal  stability/reactivaiion  kinetics 

—  H/D  in-diffusion  SIMS  elemental  depth  profiles,  effects  of  variously  doped  epi-over  layers 

—  correlation  of  elemental  and  electrically  active  depth  profiles,  influence  of  type  of  plasma 
exposure 

Passivation  of  deep  levels: 

—  passivation  of  electron  trap  deep  levels  in  MBE  GaAs,  shown  by  DLTS 

—  thermal  stability  exceeding  that  of  shallow  level  passivation 

—  associated  photoluminescence  efficiency  increase,  particularly  in  AlGaAs  capped  layers 
AIGaAs 

Passivation  of  donors: 

—  passivation  of  Si  and  Te  donors 

—  thermal  stability/reactivaiion  kinetics  of  Si  donors 
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Passivation  of  deep  levels: 

—  passivation  of  the  DX  center 

—  thermal  stability/rcactivation  kinetics  of  DX.  shown  to  be  basically  the  same  as  for  Si 
donors  except  for  exhibiting  a  distribution  of  dissociation  energies 

—  supported  model  of  an  isolated  donor  in  a  non-unique  site,  in  a  region  of  inhomogeneous 
alloy  composition 

GaAs-on-Si 
Passivation  of  defects: 

—  increase  in  the  reverse  breakdown  voltage  of  Schoniy  diodes 

—  only  partial  passivation  of  interface -originating  dislocations 

—  comparison  of  the  effectiveness  of  H  diffusion  and  H  implantation  +  anneal 

InP 

Passivation  of  acceptors: 

—  evaluation  of  various  H  permeable  surface  protection  layers 

—  strong  passivation  of  acceptors 

—  thermal  stability/reactivation  kinetics 

—  SIMS  depth  profiling  of  H/D  diffusion  into  p-InP  and  correlation  with  passivation  depth 

—  phoioluminescence  study  of  passivated  p-InP 

ELEMENTAL  semiconductors 
Si. 

Shallow  level  passivation: 

—  passivation  of  thermal  donors 

—  vibrational  spectroscopy  of  acceptor-H  and  donor-H  complexes  in  Si 

—  thermal  stability  of  donor-H  complexes 

—  uniaxial  stress  determination  of  the  symmetry  of  donor-H  and  acceptor-H  complexes  in  Si. 

—  determination  of  the  kinetics  of  the  motion  of  H  between  bond  -  centered  sites  around  the 
B  atom  in  the  B-H  complex  in  Si 


HYDROGEN  DIFFUSION  AND  DOPANT  NEUTRALIZATION  IN  1 1 1 - V  CONFOUNDS . 


J.  CHEVALLIER,  B.  THEYS ,  A.  JALIL,  R.  RAHBI ,  J-C.  FESANT 
Laboratoire  de  Physique  des  Solides  de  Bellevue 
1  place  Aristide  Br  iand,  92195  MEODON  (France) 


TOPICS  OF  INTEREST: 

.  We  are  currently  studying  the  hydrogen  diffusion  process  and  the  neutraliza¬ 
tion  of  shallow  dopants  by  hydrogen  in  various  III-V  crystalline  compounds 
and  alloys  (GaAs,  InP,  GaAlAs). 

These  studies  are  very  usefully  completed  by  high  resolution  infrared  opti¬ 
cal  spectroscopy  to  get  information  on  the  microscopic  description  of  the  H- 
dopant  complexes  (Coll,  with  B.PAJOT,  GPS-ENS  Paris). 

.  Applications  of  the  dopant  neutralization  by  hydrogen  are  investigated  in 
collaboration  with  other  laboratories  (for  the  HFET  see  the  work  presented  by 
E.  CONSTANT  and  N.  CAGLIO  from  the  University  of  Villeneuve  d'Ascq). 

MATERIALS- 

-  LPE  grown  GaAs:Si(p+)  (M.C.B01SSY,  BTC  Compelec,  CAEN) 

-  MOCVD  and  MBE  grown  GaAsiSi  (n+)  and  GaAlAs : Si  (n)  ( R .  AZOULAY ,  CNET  BA- 
GNEOX;  C.W.TO,  ATT  Bell  Laboratories,  MUBRAY  HILL) 

-  MOCVD  grown  GalnAs/InP:  Zn  structures  (B.  ROSE,  CNET  BAGNEDX) 

HYDROGEN  T)I FFIIST ON  IN  GaAs  AND.  InP;. 

.  Hydrogen  is  introduced  from  a  R.F.  hydrogen  plasma: 

Ts a e p l e  :  100°C  -  300°C ,  R.F.  power  -  0.01-0.04  W/cm3 
.  Hydrogen  detection  is  performed  by  SIMS  analysis  (CAMECA  IMS  3f  system). 

E&s.ults ;  The  close  correspondence  between  the  net  acceptor  concentration  and 
the  hydrogen  concentration  in  the  plateau  region  for  diffusion  temperatures 
up  to  300“C  [1,2]  is  a  signature  of  a  H-acceptor  pair  formation  (Fig.  1  and 

2).  Its  existence  up  to  300“C  supports  the  idea  of  a  relative  thermal  stabili¬ 
ty  of  these  complexes  and  then  suggests  that  H~acceptor  interactions  have  to 
be  taken  into  account  in  the  hydrogen  diffusion  modelling: 


Fig.  1  -  D  concentration  profiles  in  p* 
type  GaAs:  Si  ( ps 1 0 1 9cm~ 3 )  exposed  to 
R.F.  D  plasma  for  90  min. 


Fig.  2  -  D  concentration  profiles 
in  a  Ga InAs/p- InP :  Zn  structure 
exposed  to  R.F.  D  plasma  for  20  min 
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In  n-type  GaAsSi,  the  hydrogen  diffusion  study  has  been  performed  between 
200  and  iOO'C.  These  profiles  are  close  to  an  erfc  function.  The  diffusion 
coefficient  has  an  activation  energy  of  1.38eV. 

DOPANT  NFDTRAI  77.ATTQN; 

Resul  i,s  :  In  p-type  GaAs:Si  and  p-type  InP:Zn,  the  electrical  activity  cf 
the  acceptors  is  significantly  reduced  after  hydrogen  diffusion.  This  is  ac¬ 
companied  by  an  increase  of  the  free  hole  nobility  (fig, 3)  [2]. 


Fig.  3  -  Hole  nobility  vs  T  curves  for 

an  InP:Zn  epilayer  before  (a)  and  after 
(b)  hydrogen  plasma  (6h , 200”C) followed 
by  an  annealing  ( 5min . , 275*0) 
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Fig.  4  -  Temperature  dependence  of 
n  and  vn  in  a  Gao . 7 s Alo . 2 sAs : Si 
alloy  before  and  after  hydrogen 
plasma  exposure. 


The  transport  properties  of  GaAs:Si  (n)  and  GaAlAs:Si  (a)  after  hydrogena¬ 
tion  show  a  drastic  reduction  of  the  shallow  donor  concentration.  Moreover, 
in  the  GaAlAs.Si  alloys,  D-X  centers  are  also  neutralized  (fig. 4)  in  agree¬ 
ment  with  previous  DLTS  experiments  [3,4].  These  D-X  centers  are  more  effi¬ 
ciently  neutralized  than  the  shallow  donors  for  x=  0.25  so  that,  in  such  a 
passivated  alloy,  the  transport  properties  are  mainly  governed  by  residual 
shallow  donors  [4],  An  increase  of  the  electron  mobility  is  observed  after  hy¬ 
drogenation  as  in  GaAs.Si(n). 

The  increase  of  the  free  carrier  mobility  in  p-type  and  n-type  materials  is 
due  to  a  reduction  of  the  ionized  impurity  concentration  due  to  the  transfor¬ 
mation  of  dopants  into  neutral  H-impurity  complexes.  In  hydrogenated  InP:Zn, 
the  hole  mobility  is  somewhat  smaller  than  the  hole  mobility  in  hydrogen  free 
material  with  similar  net  acceptor  concentrations.  The  hole  mobility  limi¬ 
ting  factor  in  InP:Zn,H  is  partly  attributed  to  the  concentration  of  ionized 
donors  present  in  the  starting  material.  This  assumption  i3  supported  by  the 
absence  of  hydrogen  passivation  of  ionized  silicon  donors  in  p-type  GaAs:Si 
[1]. 
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CHARGE  STATES  CF  MIGRATING  HYDROGEN  IN  SEMICONDUCTORS*’ 

N.  M.  Johnson^’ 

Institute  of  Applied  Physics 
University  of  Erlangen-Nurnberg 
D-8520  Erlangen,  F.R.G. 

A  fundamental  issue  that  continues  to  challenge  our  under¬ 
standing  of  hydrogen  in  crystalline  semiconductors  is  the  question 
of  whether  migrating  hydrogen  can  exist  in  different  charge  states. 
Given  below  is  a  summary  of  the  experimental  evidence  that  bears  on 
this  issue,  with  emphasis  on  the  most  recent  studies;  calculations 
of  the  electronic  states  of  H  in  semiconductors  are  presented  else¬ 
where  in  this  workshop.  (Here,  the  symbols  H~,  H°,  and  H-  designate 
charge  states  of  possible  complexes  involving  interstial  H  and  host 
atoms.)  To  date,  the  subject  has  been  addressed  almost  exclusively 
in  Si.  However,  t.ie  techniques  and  issues  identified  in  Si  should 
be  applicable  to  compound  semiconductors . 

Depth  profiles  of  diffused  H  or  D  in  Si  display  several 
distinctive  features.  In  p-type  Si  the  D  profile  typically  displays 
a  plateau  at  the  uniform  B  concentration  with  the  depth  of  D 
penetration  equaling  the  depth  of  B  neutralization  [1],  No  such 
plateau  accompanies  shallow-donor  neutralization  in  n-type  Si  [2], 
The  D  closely  tracks  a  nonuniform  B  profile  in  p-type  Si  but  not  in 
counter-doped  n-type  Si  [1].  Such  features  have  been  interpreted 
with  the  hypothesis  that  H  possesses  a  deep-donor  level  within  the 
bandgap  so  that  H  migration  in  p-type  Si  includes  electric- 
field-assisted  diffusion  of  H-  and  neutralization  of  boron  proceeds 
by  the  reaction  H*+B-  ->  (HB)°  [3,4]. 

The  dramatic  effect  of  a  space-charge  layer  on  B  neutralization 
in  Si  has  also  been  interpreted  in  terms  of  the  deep-donor  model 
[3-6].  Both  in  Schottky  diodes  that  were  bias-temperature  stressed 
after  hydrogenation  [5]  and  in  n~-p  junction  diodes  that  were 
reverse-biased  during  hydrogenation  [6]  it  was  shown  that  B 
neutralization  occurs  preferentially  in  the  p-region  at  the  edge  of 
the  bias-dependent  depletion  layer.  While  this  can  be  ascribed  to 
the  drift  of  H~  in  the  depletion  layer,  it  does  not  by  itself 
establish  the  existence  of  mobile  H*  since  B  neutralization  could, 
in  principle,  involve  H°  and  free  holes  [6]. 

Studies  of  hydrogen  immobilization  in  Si  p-n  junctions  have 
provided  further  evidence  for  the  deep-donor  model  [7-9].  The 
space-charge  layer  was  used  to  reveal  charge-state  or  Fermi-level 
effects  during  hydrogen  migration  by  variations  in  the  depth 
profile  of  the  total  H  concentration.  In  n~-p  junctions,  the 
distribution  is  highly  structured  with  features  that  depend  on  a 
reverse  bias  applied  during  H  migration.  The  features  represent 
accumulations  of  H  in  excess  of  the  B  concentration  and  must 
involve  the  migration  of  an  H*  species .  Most  of  the  accumulation  is 
in  the  form  of  highly  immobile  neutral  entities,  which  are  proposed 
to  be  hydrogen  pairs,  Ha.  Formation  of  Ha  by  the  conventional 
reaction  2H°  ->  Ha  is,  by  itself,  incapable  of  accounting  for  the 

structure  in  the  depth  profiles.  However,  the  results  can  be 


explained  by  postulating  that  the  competing  reaction  H~-H°  -> 

Ha-h-,  where  h*' denotes  a  free  hole,  becomes  dominant  in  p-type  Si. 

Depth  profiles  of  D  in  uniformly-doped  n-type  Si  and  in 
epitaxial  n-n  layers  suggest  that  H  may  also  possess  a 
ceep-acceptor  level  within  the  £i  bandgap  [S]. 

In  a  new  approach  to  the  study  of  hydrogen  migration, 
electrical  measurements  are  now  being  conducted  during 
hydrogenation  [8],  The  principal  value  of  such  "real-time" 
measurements  is  the  ability  to  study  the  charged  species  as  it 
migrates  through  the  material  rather  than  attempting  to  infer  the 
charge  state  of  the  migrating  species  from  the  properties  or 
formation  kinetics  of  the  various  stable  H-related  complexes.  To 
illustrate,  the  resistivity  of  p-type  Si  (SOS)  was  observed  to 
increase  with  time  after  exposure  to  monatomic  hydrogen  at  300CC. 
If  the  ( HB ) °  complexes  are  completely  dissociated  at  300°C,  the 
increase  in  resistivity  should  reflect  compensation  due  to 
in-diffusion  of  H~ .  Similar  results  have  been  obtained  on  n-type  Si 
(SOS)  at  300°C  [10],  which  suggests  compensation  by  migrating  H~  if 
H-donor  complexing  can  be  neglected.  Clearly,  in-situ  measurements 
can  provide  direct  access  to  the  properties  of  migrating  hydrogen. 

In  compound  semiconductors  the  plateau  in  the  depth  profile  of 
diffused  D  at  the  uniform  acceptor  concentration,  which  has  been 
observed  in  p-type  GaAs : Zn  [11],  p-type  GaAstSi  [12],  and  p-type 
InP:2n  [13],  provides  the  only  available  evidence  to  suggest  a 
bandgap  donor  level  for  H  in  these  materials.  It  is  anticipated 
that  this  situation  will  soon  change. 

*>  Work  supported  by  the  Alexander  von  Humboldt  Foundation,  F.R.G. 
Permanent  address:  Xerox  PARC,  Palo  Alto,  CA  94304,  U.S.A. 
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IMPURITY  AND  STRUCTURAL  STUDIES  OF  HYDROGENATED  GaAs  AND  Ga.As.Si 
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Department  of  Electrical  and  Computer  Engineering 
I’nivcristy  of  Illinois  at  t'rbar.a-Ch.impaign,  Urbana,  IL  61801 

Hydrogenation  of  GaAs  and  AlOaAs  is  found  to  result  :n  a  significant  change  of  the  electrical  And 
optica!  properties.  Much  work  has  been  focused  on  die  hydrogenation  of  various  dopants  in  Ga.As  and 
AlOaAs.  How  ever,  there  have  been  few  reports  or.  the  scroctuial  and  optical  properties  of  these  materials 
after  hydrogenation. 

In  this  work,  the  effects  of  hydrogenation  on  the  residual  impurities  in  high  purity  n  and  p-type 
GaAs  win  be  reviewed.  The  samples  were  chnrarterirrri  before  and  after  hydrogenation  using 
complementary  characterization  techniques;  phototherma!  ionization  spectroscopy  (PTiS).  tnw  ler-pernture 
photolumincjccncc  (TL).  capacitance-voltage  (C-V).  ami  Hall-effect  measurements.  The  effects  ot 
hydrogenation  and  subsequent  annealing  on  GaAs  grown  on  Si  have  been  characterized  by  C  V,  Hall-effect 
measurements,  and  transmission  electron  microscopy  (TEM).  Previous  results  in  which  hydrogenation 
was  used  to  realise  isolation  regions  in  single  and  multiple  stripe  geometry  lasers  will  also  be  uesc.-ibed. 

After  hydrogenation,  die  concentration  of  Si  donors  and  Si  acceptors  that  are  present  in  GaAs  show 
a  decrease.  Carbon  acceptors  in  both  GaAs  and  AlGaAs  show  a  large  decrease.  TEM  reveals  that  the 
hydrogen  plasma  slightly  etches  die  surface  of  GaAs  3nd  energy  dispersive  speetocscopy  shows  an  arsenic 
deficient  surface.  The  effects  of  hydrogenation  is  observed  io  be  more  dramatic  on  GaAs  grown  on  Si. 

It  is  speculated  that  twins  and  dislocations  are  responsible  for  providing  a  path  for  significant  hydrogen 
accumulation.  The  lasers  produced  using  the  hydrogenation  process  are  uniform  with  low  threshold  currents 
and  cw’  room  temperature  operation. 
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Chaicogen  double  donors  in  Si 

Telleriura,  selenium,  and  sulfur  form  isolated  defects  and 
pure  impurity  pairs  in  n-type  silicon  which  act  as  double 
donors.  Their  activation  energies  have  been  determined  by 
Kail  effect  - , DLTS  and  IR  absorbtion  measurements  /I/. 

The  electrical  activity  of  all  these  double  donors  can  be 
completely  removed  by  hydrogen  passivation  /2,3/.  The  hydro¬ 
genation  was  performed  by  a  remote  microwave  plasma  e.g. 
at  1 50°C  for  30  min.  The  reactivation  of  the  donor  activity 
was  studied  for  the  neutralized  sulfur  double  donors  by 
the  rapid  isothermal  annealing  technique.  The  reactivation 
energies  for  the  isolated  S-impurity  and  for  the  pure  S-pair 
have  approx,  the  same  values;  they  result  in 

E.  (S°/+)  =  E  (S°/+)  =  2, it  0,15  eV  /3/. 

d  a  / 

From  the  annealing  kinetics,  it  appears  that  a  single  hydrogen 
atom  bonded  at  the  defect  site  is  sufficient  to  remove  both 
energy  levels  of  the  double  doner  center.  The  large  reactivation 
energies  suggest  that  the  hydrogen  is  strongly  bonded  at  the 
defect  site. 

Yapsir  and  coworkers  /4 /  have  recently  performed  cluster 
calculations  of  isolated  sulf ur-hydrogen-complexes  in  silicon. 
They  predict  two  stable  configurations  containing  one  and 
two  hydrogen  atoms,  respectively.  The  hydrogen  atoms  bond  to 
the  nearest  silicon  atoms  and  interact  only  weakly  with  the 
sulfur.  Already  the  inclusion  of  one  hydrogen  atom  moves  the 
ground  levels  of  the  double  doner  complex  below  the  valence 
band  edge.  DLTS  investigations  are  therefore  not  suited  to 
decide  which  of  the  two  stable  configurations  is  advantageously 
formed  by  the  hydrogenation  process.  Optical  investigations 
are  in  preparation  to  clarify  this  open  question. 


Cxyger.-related  cefects  in  St 

In  Cz-gr own  silicon,  the  oxygen  content  -ay  reach  ccncer.traticr.s 
oetween  5*10*^  and  2-10'®  cm  ,  Tr.e  cxycen  is  mainly  incorporated 
into  the  Si  lattice  cn  electrically  inactive  interstitial  sites. 
During  heat  treatments,  cxygen  ferns  clusters  or  precipitates 
which  act  as  donor  traps.  "Thermal  Donors"  (TD'c)  and  "New  Cxycer. 
Donors"  (ND's)  are  generated  in  the  temperature  ranee  of  350  - 
500°C  and  550  -  800°C,  respectively.  The  TD’s  appear  as  discrete 
peak  in  the  DLTS  spectrum  at  approx.  T  =  55  K.  Hydrogenation 
at  120  -  30Q°C  for  30  min  leads  to  a  maximum  reduction  cf  the 
TD  concentration  by  a  factor  of  approx.  5  / 5/.  This  amount  of 
TD  passivation  is  not  a  drastic  effect  and  does  not  support 
those  TD  models  which  explain  the  TD  activity  by  broken  Si  bonds. 

The  DLTS  spectra  of  ND  traps  indicate  the  presence  of  deep 
levels  that  are  continuously  distributed  in  energy.  Intensive 
DLTS  studies  have  been  performed  on  the  passivation  cf  ND  traps 
/ 5 , 6 / .  it  has  been  demonstrated  that  energy  levels  with 
Ec  -  E  >,  200  meV  (  Ec  =  conduction  band  edge  )  can  be  completely 
removed  and  that  energy  levels  with  energies  Ec  -  E  4  200  neV 
are  only  partially  reduced.  This  oehavior  supports  our  suggestion 
that  the  continuous  energy  distribution  of  the  ND  traps  is 
composed  of  two  different  types  cf  energy  states,  we  suggest 
that  interface  states  at  the  surface  of  SiOx-precipitates  as  well  as 
bound  states  in  the  Coulombic  wells  of  fixed  positive  charges 
which  reside  in  the  SiO^-precipitates  contribute  to  the  ND  traps. 

References  : 

1.  G. Rensl , G. Roos , C . Holm, F .Wagner  :  Defects  in  Semicond . , ed . 

H.J.von  Bardeleben,  Materials  Science  Forum  10-12  |1986)p.911 

2.  G . pensl , G . Roos , C . Holm, E. Sirtl , N .M .Johnson  :  Appl. Phys. 

Lett.  51 (1 987) ,451 

3.  G. Pensl , G. Roos , P . Stolz , N . M . Johnson ,C . Holm  :  Mat.Res.Soc. 

Symp.Proc. 104(1988) ,241 

4 .  A. S. Yapsir , P. Deak, R.K. Singh, L.C. Snyder, J.W. Corbett, T.M.Lu  : 
submitted  to  Phys.Rev.B,  April  1988 

5.  K .  Holzlein ,  G .  Pensl  ,M .  Schulz  ,N  Johnson  :  Appl .  Phys  .  Lett . 

48  (  1  986  )  ,  9  1  6 

6  .  G . Pensl  ,i4 . Schulz , K .  Holzlein , w.  Serghol z ,  J . L.  Hutchison  : 

Appl. Phys. A,  Nov. 1988 


_  30- 


HYDROGEN  PASSIVATION  OF  DOPANTS  AND 
DEFECTS  IN  lll-V  COMPOUNDS  AND  THEIR  ALLOYS 


Extended  and  point  defects 


J.  H.  Neethling  and  H.  C.  Snyman  (University  of  Port  Elizabeth):  "TEM  study  of  the 
extended  crystal  defects  in  hydrogen-implanted  GaAs ".  Invited  talk. 

P.  Hautojarvi  (Helsinki  University  of  Technology)  and  C.  Corbel  (C.E.N.  Saclay): 

Positron  annihilation  spectroscopy  of  defects  in  as-grown  and  irradiated 
GaAs".  Invited  talk. 

J.  M.  Zavada  (E.R.O.  London),  H.  A.  Jenkinson  (US  Army  ARDEC  Dover),  R.  G. 
Wilson  (Hughes  Malibu)  and  S.  J.  Pearton  (A.T.&T.  Bell  Labs  Murray  Hill):  "Dopant 
neutralization  and  defect  decoration  in  hydrogen-bombarded  gallium 
arsenide".  Invited  talk. 
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TEM  STUDY  OF  THE  EXTENDED  CRYSTAL  DEFECTS  IX  HYDROGEN 
IMPLANTED  GaAs 


J  H  NEETKLING  and  H  C  SKYMAN 


Department  cf  Physics,  University  cf  Port  Elizabeth, 

P  0  Sex  1600,  Port  Elizabeth,  6000,  South  Africa 

Proton  bombardment  of  GaAs  has  become  a  standard  technique  in  the 
fabrication  of  various  optoelectronic  devices.  Because  the  point 
and  extended  defects  introduced  curing  bombardment  critically 
influence  the  electrical  properties  of  the  device,  the  knowledge 
of  the  nature  and  distribution  of  the  implantation  damage  is 
important . 


Cross-sectional 

transmission 

electron 

microscopy  was  used 

to 

investigate  the 

crystal  structure  at 

and  near  the  end  cf 

the 

projected  range 

cf  300  keV 

protons 

implanted  at  T°  eff 

the 

GaAs  [OCl]  to  total  doses  cf  5xl015,  lxlQ16  and  lxlO1'  H  cm  2. 

In  the  low  dose  (1015-1016  H+  cm  2)  samples,  no  visible  defects 
were  observed  without  annealing.  After  a  15  minute  anneal  at 
SOO^C,  defects  in  the  form  of  black  dots  became  visible.  With 
further  increase  in  temperature  (500  to  vi  (  rw<-.r_.  •=-<-- 

grew  in  size  and  were  analysed  to  be  hydrogen  platelets  (i.e. 
hydrogen  filled  vacancy  loops)  on  the  (110)  planes  cf  GaAs.1'2 

15  +  -  2 

In  the  5x10  H  cm  samples  annealed  at  800  °C  and  abeve  cr  the 
16  +  —  2 

10  H  cm  samples  annealed  at  55  °C  ar.d  above,  the  damage 
precipitation  was  heterogeneous  consisting  not  only  cf  platelets 
but  also  containing  damage  rafts  (i.e.  three  dimensional 
dislocation  -  void  clusters)  .  1  The  damage  rafts  were  found  to 
consist  of  planar  arrays  cf  voids  on  the  (110)  planes  and  were 
elongated  along  the  < 1 1 0 >  directions.  It  was  proposed  that  damage 
rafts  originate  at  microsplits  cn  the  (110)  cleavage  planes  of 
GaAs  following  the  cracking  open  cf  small  hydrogen  filled 
platelets  when  the  internal  gas  pressure  exceeds  that  which  is 
necessary  for  crack  propagation. 


By¬ 


line  damage  rafts  were  fcur.d  to  act  as  dislocation  sources  with  the 
generation  of  glissile  dislocations,  cf  the  type  (a/C)  <liC>  or. 
the  (111}  slip  planes  intersecting  the  rafts  cr.  the  {11C}  planes, 
in  a  similar  manner  as  usually  demonstrated  for  a  Frank-Read 
source.  The  clissile  dislocations  generated  at  the  raffs  were 
found  to  glide  over  distances  of  about  1i.it,  during  the  15  minutes 
cf  anneal  at  temperatures  in  the  range  of  500  to  800  CC. 

In  the  high  dose  samples  (10^  K  cm  2)  the  radiation  damage  was 
detectable  in  the  T  EM  without  annealing.  Upon  annealing  the 
diffuse  damage  structure  changed  quite  markedly  with  increasing 
temperature  and  at  580  °C  hydrogen  platelets,  voids,  line 

4 

dislocations  and  extrinsic  Frank  loops  became  visible. 

In  conclusion  it  is  suggested  that  upon  implantation  the  hydrogen 
atoms  diffuse  inter stitially  until  trapped  at  lattice  vacancies 
created  during  implantation.  Upon  subsequent  annealing,  at 
temperatures  of  200  °C  and  above,  the  hydrogen  atoms  become  mobile 
and  leave  the  vacancies.  At  500  “C,  the  temperature  where  vacancy 
coalescence  starts  to  take  place,  the  platelets  are  formed  by  the 
cooperative  coalescence  of  vacancies  and  hydrogen  atoms.2 

.Models  will  be  presented  to  explain  the  absorption  of  hydrogen 
atoms  at  platelets  and  to  calculate  gas  pressures  in  platelets. 
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Positron  Annihilation  Spectroscopy  of  Defects  In  As-Grown  and 
Irradiated  GaAs 
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Positron  lifetime  spectroscopy  give  information  on  vacancy-tyoe 
defects  in  the  bulk.  No  positron  trapping  is  observed  in  as-grown  Si, 
whereas  In  lll-V  compounds  positrons  have  revealed  high 
concentrations  (>1016  cm*3)  of  native  vacancy  defects.  In  LEC-grown 
GaAs  positron  trapping  in  monovacancy  defects  is  seen  in  n-type,  but 
not  In  p-type  nor  In  semi-insulating  crystals.  The  defects  are  removed 
by  annealing  above  700  K.  They  are  identified  as  As  vacancies.  Two 
Fermi-level  controlled  transitions  have  been  found  in  positron  lifetime 
spectra  at  0.032  eV  and  0.10  eV  below  the  conduction  band.  They  have 
been  attributed  to  (2-1  -)  and  (-/C)  ionisation  levels  of  the  As  vacancy. 

A  beam  of  low-energy  positrons  can  be  used  to  profile 
vacancy-type  defects  in  near-surface  regions.  This  technique 
compmed  with  ion  beam  methods  (NRB  and  RBS)  has  been  applied  to 
characterize  damage  and  hydrogen  profiles  in  n-type  Si(100)  implanted 
with  101®  35-1C0  keV  H+  ions/cm^.  The  damage  distributions  do  not 
correlate  with  the  deposited  energy  distribution.  A  vacancy-type 
defected  ever, 'aye:  is  seen  close  to  the  surface  well  above  the  range  cf 
H  atoms.  Similar  experiments  are  going  on  in  hydrogen-implanted 
semi-insulating  GaAs(100)  and  preliminary  results  will  be  discussed 
;n  the  workshop. 
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Recently,  considerable  interest  has  focused  on  the  effects  cf  hydrogen 
(H)  in  C-aAs  since  it  has  been  established  that  H  can  passivate  various 
shallow  and  deep  impurity  dopants  in  the  crystal.  Experiments  have  been 
performed  to  introduce  H  into  device  structures  to  improve  electronic 
performance  and  into  crystal  growth  procedures  to  produce  higher  quality 
GaAs  epilayers.  Implantation  of  H  into  r.-type  or  p-type  GaAs  is  a  related 
processing  technique  that  has  been  useful  for  electrical  device  isolation , 
surface  passivation,  and  optical  index  definition.  However,  most  implant¬ 
ation  studies  have  centered  on  damage  defects  induced  by  the  implantation 
and  not  on  the  additional  effects  due  to  the  chemical  activity  cf  H  in  the 
crystal. 

In  a  trier  study,  we  have  examined  the  pcst-impiar.tation  behavior  cf  H 
in  r.+  GaAs  (Si)  using  secondary  icn  mass  spectrometry  (SIMS)  and  infrared 
reflectance  measurements.  It  was  shewn  that  implanted  H  is  redistributed 
with  post-implantation  furnace  annealing  and  produces  optical  effects 
associated  with  carrier  removal  (11. 

Here  we  first  report  on  the  properties  of  H  implanted  into  p+  GaAs(Zn) 
crystals  and  the  subsequent  behavior  with  post-implantation  annealing. 
Depth  profiles  for  both  H  end  Zn  atoms,  as  a  function  cf  anneal  temper¬ 
ature,  have  been  determined  using  SIMS.  Infrared  reflectance  measurements 
indicate  little,  if  any,  change  neutralization  associated  with  H  migration 
into  the  substrate.  Results  of  these  experiments  with  p+  GaAs(Zn)  are 
compared  with  data  obtained  for  H  in  nm  GaAs  (Si)  and  with  related  studies 
involving  plasma  hydrogenation. 

The  second  part  of  our  presentation  concerns  the  redistribution  cf 
implanted  H  ions  in  GaAs  multilayer  structures  due  to  pcst-implar.taticn 
annealing.  Several  different  multilayer  structures  have  been  investigated 
including  GaAs-on-Si  samples  and  GaAs/AlAs  superlattices.  In  each  experi¬ 
ment,  SIMS  measurements  show  that  post-implantation  annealing  causes  the  H 
atoms  to  undergo  a  redistribution  and  to  accumulate  at  interface  regions 
within  the  crystal. 

For  the  GaAs-cn-Si  samples,  a  major  accumulation  cf  migrating  hydrogen 
has  been  measured  at  the  GaAs/Si  hetero-^r.terfye  J2).  This  accumulation 
reaches  concentrations  in  excess  cf  10  H  cm  and  is  stable  even  at  a 
temperature  of  700°C.  Kith  materials  having  such  a  large  lattice  mismatch 
(<*),  it  can  be  expe.ted  that  the  interface  region  would  contain  a  high 
density  of  defects  and  dangling  bends  leading  to  trapping  cf  H  atoms. 
However,  the  experiments  with  the  GaAs/AlAs  superlattices  shew  that  even  a 
hcmc-expitaxial  interface  can  produce  similar  effects  [3).  While  the 


effects  [3]  .  While  the  numerous  material  interfaces  in  these  multilayer 
structures  affect  the  migration  cf  the  atoms,  the  GaAs  buffer  iayer/C-aAs 
substrate  interface  appears  to  be  the  most  important.  This  interface 
impedes  the  migration  cf  ^pms  ir,to  the  substrate  ar.d  retains  them  at 
concentrations  cf  about  10  H  cm  ' ,  even  after  annealing  at  50C°C.  lr, 
either  experiment,  the  accumulation  of  H  appears  to  be  due  to  the  bcr.di.-.g 
cf  H  with  dislocations  and  crystal  imperfections  created  during  the  initial 
stages  of  epitaxy. 

These  findings  are  used  to  form  a  qualitative  model  for  the  behavior 
of  H  implanted  into  GaAs.  Initially,  the  majority  cf  H  decorates  crystal 
damage  within  the  implanted  region.  During  pcst-implanaticn  annealing,  H 
atoms  are  released  from  the  damage  centers,  migrate  to  other  regions  of  the 
crystal,  and  are  trapped  either  by  different  damage  centers  or  by  impurity 
atoms.  Charge  neutralization  is  due  to  ion  irradiation  effects,  as  well  as 
to  trapping  of  H  at  impurity  dopant  locations.  The  redistribution  cf 
implanted  H  is  a  complicated  process  depending  on  many  parameters  including 
damage  concentration,  dopant  level  and  dopant  species.  Further  experiments 
are  needed  to  clarify  the  relative  importance  cf  these  parameters. 


II]  L.L.  Liou  et  al,  "Infrared  Study  of  Hydrogen-  and  Carbon-implanted 
Heavily  Doped,  n-Type  GaAs",  J.  Appl.  Phys.  5%  1936  (1986). 

[2]  S.J.  Pearton  et  al,  "Effects  of  Atomic  Hydrogen  Incorporation  in 
GaAs-on-Si",  J.  Appl.  Fhys.  (To  be  published). 

[3]  J.M.  Zavada  et  al,  "Redistribution  of  H  and  Se  in  GaAs/AlAs  Multilayer 
Structures  with  Post-Implantation  Annealing",  J.  Appl.  Phys.  (To  be 
published) . 
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B.  Cletjaud,  D.  Cote,  M.  Krause,  C.  Porte  (Universite  Pierre  et  Marie  Curie  Paris)  and 
W.  Ulrici  (Zentralinstitut  fur  Elektronenphysik  Berlin):  "Spectroscopic  studies  of 
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B.  Pajot  (Universite  Paris  7):  "Spectroscopic  investigations  of  the 
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J.  Tatarkiewicz  and  M.  Stutzmann  (MPI-FKF  Stuttgart):  "Hydrogen  vibrations  in 
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The  lattice  locations  of  isolated  substitutional  dopant  atoms  with  a  low 
mass  in  GaAs  can  be  determined  by  measurements  of  their  infrared  localized 
vibrational  mode  (LVM)  absorption  using  high  resolution  Fourier  transform 
spectroscopy.  It  is  a  good  approximation  to  assume  that  second  nearest 
neighbours  are  stationary  so  that  the  angular  frequency  of  a  mode  may  be 
written  as  u2=k[Mlmp*(xMnn)"1  ] ,  where  Mimp  is  the  mass  of  the  impurity,  Mnn  is 
the  mass  of  the  nearest  neighbours  and  the  parameter  x  depends  on  the  local 
bending  and  stretching  force  constants  There  is  a^arge  impurity  isotope 
shift,  (e.g.  between  say  2SSi  and  Si  and  'h  and  D  but  there  Eire  also 
measurable  frequency  changes  when  nearest  neighbour  isotopes  are  changed. 
Atoms  occupying  Ga  lattice  sites  (As  neighbours  :  As,  100%  abundant)  produce 
sharp  lines,  whereas  broadened  lines  with  structure  occur  for  atoms  occupying 
As  lattice  sites  (Ga  neighbours:  6  Ga,  60%  abundant  and  7IGa  40%  abundant). 

When  impurity  hydrogen  pairs  form  the  symmetry  is  reduced  from  Ta  (to  say 
Cdv)  so  that  the  impurity  will  have  a  longitudinal  mode  along  the  pair  axis 
and  a  doubly  degenerate  transverse  mode  ui.  The  modified  local  bonding  will 
change  the  local  force  constants  :  a  reduction  in  the  frequency  of  an  LVM 
indicates  a  corresponding  reduction  in  the  force  constants.  It  should  be 
remembered  that  stretching  force  constants  are  about  t^n  times  greater,  and  are 
therefore  more  important,  than  the  bending  constants  1  .  If  the  frequency  is 
lowered  too  much,  the  mode  will  become  incorporated  into  the  lattice  continuum 
and  will  not  be  detected.  The  impurity  modes  will  also  be  lowered  in 
frequency  when  hydrogen  is  replaced  by  deuterium.  In  general  the  presence  of 
a  static  charge  on  an  impurity  in  a  covalent  host  should  not  change  the 
strength  of  the  IR  integrated  absorption  appreciably  compared  with  that 
produced  by  a  neutral  impurity121.  Thus  reasonable  conservation  of  integrated 
absorption  is  expected  when  hydrogen-impurity  pairs  are  formed,  and  the 
strength  of  uj.  should  be  about  twice  that  of  un. 

For  the  paired  hydrogen,  longitudinal  and  transverse  modes  are  also 
expected.  The  ratio  r  of  a  hydrogen  frequency  uH  to  that  of  the  corresponding 
deuterium  mode  od  gives  information  about  the  mass  of  the  atom  to  which  the 
hydrogen  is  bonded.  Although  most  of  the  kinetic  energy  will  be  in  the 
H-mode,  the  impurity  can  also  move  and  transfer  energy  to  other  neighbours. 
It  follows  that  the  value  of  r  will  be  smaller  than  the  value  calculated  for  a 
diatomic  molecule.  Finally,  the  application  of  a  uni-axial  stress  along 
certain  directions  will  remove  the  double  degeneracy  of  both  the  impurity  and 
the  hydrogen  transverse  modes. 

These  topics  are  illustrated  for  the  systems  Bec«-H131  and  Sic»-HU1  in 

GaAs. 
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MBE  GaAs  doped  with  Be  to  a  concentrat  i  on  ogf  3.7xl0'8cm  3  was  fully 
cor.pensated  by  2MeV  electron  irradiation.  The  Beta  impurities  gave  a 
symmetrical  LVM  line  at  482.4cm’1 2 3 4 5 6  with  a  half  width  &  of  0.5cm  as  expected 
for  an  impurity  on  a  Ga  lattice  site.  The  integrated  absorption  was 
characterized  by  an  apparent  charge  of  i)=l.li0.1e,  in  agreement  with  previous 
measurements  on  LPE  material.  However,  as-grown  material  containing  Bec« 
appeared  to  show  a  much  stronger  LVM  line  with  7j=1.7e.  The  line  was  shifted 
to  a  lower  energy  by  0.3cm  1  and  most  unexpectedly  the  value  of  A  increased  to 
1.0cm’1.  These  differences  have  now  been  resolved  by  recognizing  that  the 
line  shape  is  an  asymmetric  Eano  profile.  Following  hydrogen  diffusion  from  a 
plasma  the  Be  LVM  is  reduced  in  strength  and  a  new  line  appears  at  555.7cm’1 
for  Be-H  and  at  553.6cm*1  for  Be-D.  Only  one  line  is  found,  implying  that 
either  ui  or  un  has  been  lowered  into  the  continuum.  Taking  account  of  the 
Fano  effect  the  new  line  has  about  2/3  the  absorption  lost  from  the  isolated 
Be(Td).  Only  one  high  frequency  LVM  line  from  hydrogen  is  found  for  Be-H  at 
2037.1cm*1  and  Be-D  at  1471.2cm"1,  giving  r=1.38  gwhich  implies  that  the 
hydrogen  is  bonded  to  a  host  lattice  atom  and  not  a  Be  impurity  (r=1.386  for 
Zn-H).  This  is  undoubtedly  a  stretching  mode  and  so  the  transverse  mode  must 
be  soft,  i.e.  at  a  frequency  less  than  about  300cm  Recent  calculations151 
support  these  interpretations  and  will  be  discussed  in  a  separate  paper.  It 
is  also  concluded  that  the  555.7cm"  mode  should  be  identified  with  oj.  rather 
than  U/i.  Thus  the  Be  is  tri-co-ordinated,  while  the  H  atom  bonds  to  an 
adjacent  As  neighbour.  It  appears  that  re-orientation  of  the  Be-H  pair  may 
occur  as  A  for  the  Be-H  (2037cm’1)  mode  is  greater  than  that  for  Be-D. 


For  Slca-H,  there  is  again  only  one  mode  from  the  paired  Si  at  409.95 
cm’1  (and  409.45cm*1  for  Si-D).  As  there  is  evidence  for  direct  bonding  of 
the  hydrogen  to  the  Sica  atom,  it  was  thought  at  first  that  the  local  force 
constants  might  Increase  so  that  both  u#  and  ui  would  be  greater  than 
o=384cm’1  for  isolated  Sica(Td).  Recent  theory  indicates  however  that  one  of 
the  Sica  bonds  to  an  As  neighbour  breaks,  in  which  case  would  occur  at  a 
low  frequency  and  would  not  be  observed.  Both  and  u>i  are  observed  for  the 
paired  hydrogen,  and  silicon  isotope  effects  are  evident  on  uh  .  Recently  it 
has  been  shown  that  the  line  ascribed  to  ui  is  Indeed  doubly  degenerate 
indicating  that  the  defect  has  C3v  symmetry. 
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All  the  GaAs,  InP  and  GaP  materials  grown  by  the  liquid-encapsuled 
Czochralski  technique  show  series  of  infrared  absorption  lines  in  the  2000 
cm~l  range  for  GaAs  and  in  the  2200-2300  cm-1  range  for  InP  and  GaP.  These 
lines  are  due  to  local  modes  of  vibration  of  hydrogen  related  complexes, 
involving  As-H  or  P-H  stretching  modes  of  vibration  [l].  In  indium  phosphide, 
some  of  the  complexes  involve  dopant  impurities  |2|  whereas  some  others  involve 
lattice  defects  [3],  In  gallium  arsenide,  no  correlation  could  be  found  with 
any  dopant  and  probably  most  of  the  observed  lines  are  related  with  complexes 
involving  hydrogen  and  lattice  defects  ffl.  The  acceptor  impurity  related 
complexes  are  believed  [l]  to  be  of  the  bond  center  type  first  proposed  by 
Pankove  et  al.  [5]  for  the  boron-hydrogen  complex  in  silicon.  In  order  to  learn 
more  about  the  complexes  involving  lattice  defects,  we  have  performed  infrared 
measurements  under  uniaxial  stress.  The  results  of  these  experiments  together 
with  the  other  properties  of  the  center  suggest  that  the  2001  cm"*  line  observed 
in  GaAs  is  due  to  the  binding  of  an  hydrogen  atom  by  an  As  dangling  bond  in 
a  gallium  vacancy  [6],  The  same  type  of  defect  involving  an  indium  vacancy 
in  indium  phosphide  has  a  local  mode  of  vibration  at  2202.4  cm”1  [7].  The 
microscopic  structure  of  the  other  complexes  involving  lattice  defects  is  not 
fully  understood  yet. 
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At  GPS,  we  have  primarily  studied  the  infra-red  absorption  of  the  hydrogen-related  modes  in 
in-V  semiconductors  plasma-treated  or  implanted  with  hydrogen.  The  frequencies  where 
absorption  is  observed,  eventual  natural  isotope  effects  and  the  ratio  of  the  frequencies  when  !H  is 
replaced  by  2H  (the  r-factor)  provide  information  enabling  to  identify  the  atom(s)  to  which 
hydrogen  is  bonded.  Once  done,  it  is  generally  possible  to  establish  the  structure  of  the  complex 
involving  hydrogen  and  to  relate  it  to  the  neutralization  of  the  dopant  atom.  The  structure  can  be 
also  deduced  or  confirmed  by  measurement  of  the  absorption  of  polarized  radiation  by  the  sample 
under  uniaxial  stress. 

We  summarize  here  the  main  results  obtained  using  low-temperature  infra-red  absorption  with 
Fourier  transform  spectroscopy. 

1)  Silicon  donors  in  GaAs 

a)  Neutralization  bv  hydrogen  of  n-tvpe  GaAs:Si 

High-resolution  measurements  have  shown  that  H  was  bonded  to  Si  on  the  Ga  sublattice. 
From  these  results  was  proposed  the  location  of  the  hydrogen  atom  on  the  anti-bonding  site  and 
the  neutralization  of  the  electrical  activity  by  localization  of  the  donor  electron  on  the  Si-H  bond. 
Uniaxial  stress  measurements  proved  to  be  in  agreement  with  the  trigonal  symmetry  expected 
from  this  centre  in  the  GaAs  lattice. 

b)  Proton  implantation  in  n-tvpe  GaAs:Si 

The  Si-H  bending  mode  already  observed  in  H-p!asma  neutralized  material  is  also  found  in  H" 
implanted  material  at  ambient.  The  co-existence  of  this  mode  with  modes  of  H-related  defects 
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show  that  in  H'  implanted  GaAs:Si,  a  small  pan  of  the  carrier  removal  already  arises  from,  direct 
neutralization  of  dopant  by  hydrogen.  This  fraction  increases  under  annealing  up  to  a  point  where 
trr.T.er  rcm.cvd  comes  c..Iy  from  dopant  neutraitzation. 

c)  Silicon  donors  in  Ga^Al  As 

Substitution  of  A1  to  Ga  atoms  usually  lowers  the  symmetry  of  the  dose  environment  of  the 
trigonal  centre,  hence  removing  the  two-fold  degeneracy  of  the  Si-H  bending  mode.  This  and  the 
different  splittings  from  the  different  arrangements  of  the  A1  atoms  produce  many  lines  near  the 
position  of  the  single  line  in  pure  GaAs.  Analysis  of  the  structure  informs  on  the  actual  probability 
to  find  a  given  configuration  as  compared  to  the  one  derived  from  the  AI  concentration  alone. 

2)  Acceptors  in  Gallium  Arsenide  and  Indium  Phosphide 

a)  The  formation  of  hydrogen-dopant  complexes  in  p-type  III-V  complexes  is  a  two-step 
process  resulting  from  i)  compensation  of  acceptors  by  hydrogen  and  ii)  complexing  of 
hydrogen  with  acceptor  dopant  driven  by  coulombic  interaction.  Spectroscopy  shows  that  this 
pairing  is  achieved  by  covalent  bonding  with  hydrogen  in  a  bond-centred  configuration. 

b)  Silicon  in  p-type  GaAs 

The  lines  observed  in  p-type  GaAstSi  neutralized  with  hydrogen  are  ascribed  to  the  stretchine 
mode  of  hydrogen  bonded  to  Si  in  the  As  sublattice.  The  r-factor  indicates  however  some 
coupling  between  hydrogen  and  the  nearest  neighbour  Ga  atom.  No  interaction  is  found  between 
hydrogen  and  the  residual  Si  donors,  showing  that  in  GaAs,  the  formation  of  complexes  with 
dopant  atoms  is  possible  only  when  these  atoms  are  electrically  active. 

c)  Acceptors  in  the  cation  sublattice 

The  observations  for  Zn  in  GaAs  an  in  InP  as  well  as  for  Be  in  GaAs  arc  explained  by  the 
binding  of  hydrogen  to  the  group  V  atom  nearest  neighbour  of  the  acceptor.  The  interaction 
between  hydrogen  and  the  acceptor  seems  so  weak  that  it  can  be  thought  of  as  tri-coordinated. 

d)  Temperature  effects  and  line  width 

No  temperature-dependent  structure  is  found  for  hydrogen  lines  with  all  the  acceptors  when 
temperature  is  raised.  This  could  be  interpreted  as  a  consequence  of  an  "on-line"  geometry  of  the 
H  bond  along  the  <1 1 1>  direction.  The  dependence  of  the  width  and  position  of  the  lines  on 
temperature  reflects  the  coupling  between  the  complex  and  the  lattice  and  a  pan  of  it  can  be 
explained  by  steric  effects.  An  interpretation  of  the  results  is  on  the  vay. 
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By  Ff-I?.-spectrosccpy  we  studied  130— grown  3a?  doped  in  the 
melt  with  Sc,  Ti ,  V,  Cr,  Co,  or  Hi.  In  C-a?:Ci  3  sharp  lines  have 
been  found  around  2200  on-”1  (Pig.  1)  which  are  local  vibrational 
mode  (LVH)  lines  as  proved  by  their  temperature  dependence.  These 
lines  are  interpreted  as  due  to  ?-H  centres  [l]  in  analogy  to  the 
As-K  centres  investigated  in  IBC-grown  3 a As  [2,  3].  However, 

these  H-related  lines  could  be  detected  only  in  C-aP:Ti. 

In  C-a?  doped  with  Sc,  V,  Cr,  Co,  or  Hi  on  the  other  hand 
several  LVH  absorption  lines  in  the  region  1000  cm  to  1400  cm- 
were  measured  (Fig.  2)  not  detectable  in  Ga?:?i  samples.  These 
lines  have  been  assigned  to  centres  containing  oxygen  on  inter¬ 
stitial  sites  ((h)  Cl].  There  is  no  correlation  of  the  lines  with 
the  hind  of  transition  metal  doping.  The  large  number  of  lines  is 
possibly  due  to  different  centres  of  0^^  bonded  to  ?  or  Ga  and 
completed  with  typical  residual  impurities  (C3,  3„  ,  3ir  )  or  ?. 
anti-site  defects  (all  of  them  present  in  the  samples  ~  1010 
cm  ).  The  absence  of  all  of  these  lines  in  C-a?:Ti  supports  their 
assignment  to  oxygen-related  centres  as  Ti  is  known  to  be  a  very 
effective  getter  material  for  oxygen.  A  first  quantitative  argu¬ 
ment  in  favour  of  0.  is  the  detection  of  a  weak  line  at  119S.2 
—1  1  —1 
cm  accompanying  the  most  intense  line  at  1275.3  cm  .  The  ratio 

of  their  frequencies  r  =  1.0635  is  very  close  to  the  expected  r  = 
1.0507  for  tv  <?  corresponding  ^0  line,  and  their  ratio  cf  intc-r.- 
sities  (0.0025)  is  in  good  agreement  with  the  natural  abundance 
of  1 -0  (0.0020). 
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The  possibility  of  passivating  many  points  defects  ir.  a  large  variety  of  different  semiconductors  by 
introduction  of  atomic  hydrogen  has  spurred  new  interest  in  this  simple  impurity  and  its  bonding  pro¬ 
perties.  However,  no  general  approach  to  hydrogen  bonding  in  crystals  has  been  developed  so  far.  In 
our  recent  paper  [  1  ]  we  have  shown  that  a  simple  model  of  a  harmonic  oscillator  embedded  in  a  dielectric 
medium  enables  us  to  predict  hydrogen  vibrational  frequencies  in  semiconductors  as  well  as  in  ionic 
crystals.  Here  we  present  the  next  logical  step  in  the  extension  of  our  model:  instead  of  using  molecular 
hydrogen  vibrational  frequencies  comto  predict  corresponding  hydrogen  frequencies  in  crystals  coc,  we 
estimate  (cf.  e.g.  [2])  the  force  constant/:  of  the  harmonic  potential  V'=  /cc2  from  the  dissociation  energy 
D  [3]  and  the  bond  length  a  [3]  of  the  corresponding  free  diatomic  molecule:  k~D:cd.  In  the  previous 
paper  [1]  we  used  a  correction  due  to  the  depolarizing  field  produced  by  the  vibrating  dipole  in  a  cavity 
inside  the  solid.  However,  in  the  case  of  semiconductors  this  correction  is  smaller  than  the  error  of  the 
estimate.  Hence  in  the  presentmodel  we  plot  (see  also  Table  I,  where  Acu  indicates  hydrogen  band  width) 

lab-lei 


Crystal 

Bond 

D 

kcal/mol 

a 

A 

v'D/a 

VeV/A 

env1 

cm  1 

Aco 

cm'1 

Ref. 

Si 

SiH 

71.4 

1.520 

1.158 

2042 

1996 

200 

[4] 

SiC 

SiH 

71.4 

1.520 

1.158 

2042 

2090 

150 

[5] 

CH 

80.0 

1.120 

1.664 

2861 

2920 

25 

Ge 

GeH 

66.8 

1.588 

1.072 

1834 

1886 

125 

[6] 

GaAs 

GaH 

68.0 

1.663 

1.033 

1604 

1711 

100 

[6] 

AsH 

65.0 

1.534 

1.094 

2130 

1998 

100 

ZnSe 

ZnH 

19,5 

1.595 

0.577 

1608 

1315 

150 

[7] 

SeH 

73.0 

1.470 

1.210 

2400 

2150 

50 

ZnTe 

ZnH 

19.5 

1.595 

0.577 

1608 

1215 

100 

[8] 

TeH 

48.2 

1.520 

0.951 

2250 

1655 

25 

ZnS 

ZnH 

19.5 

1.595 

0.577 

1608 

1680 

100 

[9] 

SH 

81.4 

1.346 

1.396 

2200 

- 

- 

CdTe 

CdH 

16.5 

1.780 

0.475 

1337 

1205 

100 

[10] 

TeH 

48.2 

1.520 

0.951 

2250 

1730 

25 

CdS 

CdH 

16.5 

1.780 

0.475 

1337 

1285 

150 

[11] 

SH 

81.4 

1.346 

1.396 

2200 

2450 

100 

GaP 

GaH 

68.0 

1.663 

1.033 

1604 

1849 

25 

[12] 

PH 

82.0 

1.422 

1.326 

2365 

2204 

25 

[13] 

InP 

InH 

59.0 

1.838 

0.870 

1476 

- 

PH 

82.0 

1.422 

1.326 

2365 

2204 

25 

[M] 

the  observed  hydrogen  vibrational  frequencies  in  crystals  a>cas  a  function  of  the  square  root  of  f  As 
shown  in  Fig.  1  this  procedure  gives  an  excellent  fit,  except  for  the  group  Da  metals  (Zn  and  Cd).  On 
the  other  hand,  if  we  plot  hydrogen  molecular  vibrational  frequencies  <nm  vs.  vit  (Fig.  2)  we  get  a  much 
worse  fit.  We  also  notice  that  the  molecular  frequencies  in  zinc  and  cadmium  hydrides  are  much  higher 
than  they  should  be  if  calculated  from  the  simple  model  of  the  harmonic  oscillator.  We  explain  this  fact 
by  the  metallic  character  of  these  bonds,  as  opposed  to  the  covalent  character  for  other  bonds  considered 
here.  Finally,  let  us  note  that  experimental  data  for  ZnS  has  been  obtained  for  a  crystal  grown  in  hydro¬ 
gen  atmosphere  [9],  whereas  in  all  other  materials  hydrogen  was  introduced  by  implantation. 

We  conclude  that  our  simple  model  applies  to  all  covalent  semiconductors,  permitting  easy  pre¬ 
dictions  of  hydrogen  vibrational  frequencies  in  these  crystals.  This  can  be  of  some  help  for  future  studies 
of  hydrogen  in  various  semiconductors  (see  e.g.  [15]),  where  the  vibration  energies  are  still  unknown. 
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\\  e  have  studied  the  hydrogen  passivation  of  B-doped  Si  by  means  of  infrared  reflectance 
and  Raman  scattering.  The  samples  used  were  bulk  doped  with  boron  concentration 
between  I  x  1019  and  1.2  x  103°  cm-3. 

The  IR  reflectance  of  the  doped  and  passivated  samples  was  measured  in  the  region 
of  the  tree-carrier  plasma  edge.  After  passivation,  a  shift  of  the  plasma  edge  to  lower 
wavenumbers  is  observed,  corresponding  to  a  decrease  of  the  free  hole  concentration. 
The  interference  pattern  observed  in  the  IR  spectra  of  the  hydrogenated  samples  allows 
us  to  deduce  the  thickness  of  the  passivated  region.  Moreover,  the  decay  of  this  in¬ 
terference  pattern  is  directly  related  to  the  steepness  of  the  H-profile  in  the  passivated 
samples.  With  these  three  variables  (hole  concentration  in  the  surface,  depth  of  the 
passivated  region  and  steepness  of  the  H-profile)  we  can  fit  the  experimental  reflectance 
spectra  for  different  B-concentrations  and  times  of  exposure  to  hydrogen.  A  typical 
experimental  spectrum  together  w-ith  a  calculated  fitting  is  presented  in  Fig.  1.  The  cal¬ 
culated  profiles  agree  well  with  those  obtained  by  means  of  SIMS1.  So,  IR  reflectance 
provides  a  direct  and  nondestructive  way  to  obtain  depth  profiles  in  passivated  sam¬ 
ples.  The  correlation  found  between  B-concentration  and  H-depth  indicates  that  the 
H-diffusion  is  basically  controlled  by  an  impurity  trapping  mechanism.  For  a  boron 
concentration  of  1  x  1019  cm-3,  the  diffusion  coefficient  obtained  at  ISO^C  is  about 
6  x  10'13  cm3  s'1. 

Direct  information  about  the  free  hole  concentration  can  also  be  obtained  by  looking 
at  the  Raman  active  silicon  zone-center  optical  phonon,  which  in  doped-Si  appears 
shifted  and  asymmetrically  broadened  in  comparison  to  the  undoped  semiconductor, 
showing  the  so-called  Fano  profile.  The  narrowing  of  this  profile  after  hydrogenation 
can  be  used  as  a  direct  test  to  obtain  insight  into  the  remaining  hole  concentration3. 

From  a  more  microscopic  point  of  view,  the  localized  vibrations  of  boron  and  hy¬ 
drogen  in  the  Si-lattice  are  direct  probes  for  the  bonding  structure  of  these  atoms. 
Moreover,  the  H-mode  observed  at  1905  cm"1  can  be  used  to  obtain  information  about 
the  symmetry  of  the  B-H  complex  in  the  lattice.  For  this  purpose,  we  have  analyzed 
the  Raman  scattering  due  to  H-vibrations  under  uniaxial  stress  along  different  lattice 
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Fig.  1.  Infrared  reflectance  spectrum  of  a 
passivated  B-doped  Si  sample.  Points  cor¬ 
respond  to  the  experimental  data  and  the 
continuous  line  is  calculated  using  the  indi¬ 
cated  parameters. 


Fig.  2.  Off-bond-center  hydrogen 
displacements  along  '112]  (towards 
C-points)  and  along  ,112]  (towards 
C-points). 


directions.  From  our  results,  and  mainly  from  the  line  splitting  observed  for  stress 
along  a  [100]  direction,  we  find  at  100K  a  non-trigonal  symmetry  for  the  B-H  com¬ 
plexes.  This  corresponds  to  off-bond-centered  interstitial  sites  for  hydrogen,  which  can 
be  modified  under  stress.  The  stress  dependence  observed  for  the  line  intensities  indi¬ 
cates  a  stress-induced  hydrogen  redistribution  among  non-equivalent  positions  around  a 
B-impurity3.  In  Fig.  2  we  present  a  scheme  corresponding  to  the  H-displacements  along 
,112]  directions  deduced  from  the  analysis  of  the  Raman  results.  Hydrogen  can  move 
between  different  energy  minima  with  barriers  of  less  than  300  meV,  in  agreement  with 
theoretical  calculations*. 

The  discrepancy  between  our  results  and  the  trigonal  symmetry  found  by  Bergman 
e(  al.1  at  4K  may  be  due  to  the  strong  T-dependence  of  the  H-mode,  as  well  as  to 
different  internal  strain  in  the  samples  (ion  implanted  vs.  bulk  doped). 
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Advanced  nuclear  methods 
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Over  the 

last 

five  years  much  has  been 

learned  about 

muonium  impuri 

ties 

in  semiconductors  and  it 

now  appears  very 

likely  that  most  of  this  information  will  be  directly 
applicable  to  the  understanding  cf  hydrogen  impurities. 

Muonium  may  be  formed  when  positive  muons  (>i+)  are  implanted 
into  nenmetallic  samples.  They  may  be  studied  by  techniques 
closely  related  to  magnetic  resonance  and  collectively 
referred  to  as  muon  spin  rotation  or  y.S%.  Mucnium  is  a 
neutral  atom  formed  ty  an  electron  bound  to  the  p+  and  may  be 
regarded  as  a  very  light  isotope  of  hydrogen  <ny=mp/S).  In 
semiconductors  muonium  is  quite  different  from  mucnium  in  free 
space.  Two  distinct  forms  of  muonium  are  seen1,  one  with  an 
isotropic  hyperfine  interaction  ranging  from  83*  of  the  vacuum 
value  (in  diamond)  to  27*  (in  CuClt.T.his  normal  muonium  or  Mu 
has  a  hyperfine  interaction  which  is  about  50*  of  the  vacuum 
value  in  Si  and  Ge  but  about  65*  in  GaAs  and  GaP.  One  can 
visualize  this  as  a  neutral  muonium  impurity  covalently-bonded 
to  atoms  in  the  crystal  so  that  about  half  (for  Si  and  Ge)  the 
spin  density  is  on  the  muonium  and  about  half  is  on  the  host 
atoms.  The  experimental  and  theoretical  evidence  supports  the 
view  that  for  normal  muonium  the  is  centered  at  the 
tetrahedral  interstitial  site.  In  addition,  a 
highly-anisotropic  form  of  muonium  is  observed.  It  is  called 
anomalous  muonium  or  Mu*  and  has  a  trigonally-symmetr ic 
hyperfine  interaction  whose  angular  average  is  only  a  few 
percent  of  the  vacuum  value.  Recent  studies  by  muon 
level-crossing  spectroscopy,  augmented  by  resolved  nuclear 
hyperfine  structure  in  the  muon  precessicnal  frequencies,  have 
shown  that  for  anomalous  muonium  in  silicon  the  electron  spin 
density  and  the  muon  are  centered  on  a  covalent  bond”.  Similar 
observations  have  been  made  on  GaAs3A  and  GaP4  and  presumably 
in  all  cases  cf  Mu*  the  muon  is  near  the  center  cf  a  covalent 
bond . 


The  recent  repcrt  of  an  electron  paramagnetic  resonance 
( EFR )  spectrum  (AA-3)  in  hydrogen-implanted  silicon^  was  the 
first  example  cf  EPR  for  isolated  hydrogen  in  a  semiconductor 
known  to  these  authors.  Comparison  of  the  scaled  hydrogen  and 
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"u:n  hyperfine  parameters  and  those  of  the  le&jes*  ^-Ei  for 
ohe  A  A  -  3  renter  err  Mu*  in  silicon  indicate  chat  the  hyperfine 
parameters  are  the  same  tr  v.-ithin  about  the  uncertainty  of  the 
less  accurate  EFR  results  ( the  hydrogen  values  must  he 
multiplied  by  ohe  ratio  cf  the  muon  to  the  proton  magnetic 
dipole  moments,  i.e.  scaled).  We  may  therefore  conclude  that 
this  mucnium  center  has  a  corresponding  hydrogen  center  which 
ecists  with  nearly  the  same  structure,  and  this  is  cuite 
likely  true  for  most  other  mucnium  centers  as  well.  The 
principal  effect  cf  the  large  mass  difference  should  oocur  in 
the  dynamics  cf  the  two  types  of  systems.  In  that  regard, 
metastable  forms  of  mucnium  have  been  observed  in  many 
materials.  This  includes  5i ,  3e ,  diamond,  SiC,  GaAs,  3a?, 

CuCl,  and  CuEr.  For  some  of  these  systems  (Bi,  diamond,  CuCi , 
Cu3r )  there  is  evidence  of  a  thermal  conversion  of  one  form  to 
another . 

At  present  one  may  conclude  with  considerable  confidence, 
based  on  the  measurements  described  and  upon  theoretical 
studies  which  support  these  conclusions,  that  the  stable 
structure  of  neutral,  isolated,  interstitial  hydrogen  in 
silicon  is  one  in  which  the  proton  is  at  the  bond  center  and 
the  spin  density  is  primarily  on  the  two  nearest  silicons.  A 
similar  model  can  also  be  inferred  for  GaAs  and  Gap  although 
it  has  not  been  demonstrated  that  this  is  the  stable 
structure. 


1.  For  a  general  review  through  1986  which  includes  oil  forms 
of  muonium  seen  (no  reports  of  new  mucnium  centers  have 
occurred  through  Aucust  1988)  see:  B.  D.  Patterson,  Rev. 
Mod.  Phys.  60,  69(1998). 

2.  R.  F.  Kiefl,  M.  Celio,  T.  L.  Estle,  S.  R.  Kreitzman, 

G.  M.  Luke,  T.  M.  Riseman,  and  E.  j.  Ansaldo,  Phys.  Rev. 
Letters  60,  224(1988). 

3.  R.  F.  Kiefl,  M.  Celio,  T.  L.  Estle,  G.  M.  Luke, 

S.  R.  Kreitzman,  J.  H.  Erewer,  D.  R.  Noakes,  E.  J.  Ansaldo, 
and  K.  Nishiyama,  Phys.  Rev.  Letters  Z9,  1780(1997). 

4.  R.  F.  Kiefl,  Hyperfine  Interactions  12.,  707(1986). 

5.  Yu.  V.  Gorelkinskii  and  N.  N.  Nevinnyi,  Sov .  Tech.  Phys. 
Lett.  13,  4 S ( 1 9 8 7 ) ;  V.  A.  Gordeev,  Yu.  V.  Gorelkinskii, 

R.  F.  Kcnopleva,  N.  N.  Nevinnyi,  Yu.  V.  Obukhov,  and 

V.  G.  Firsov,  preprint  1340  from  Leningrad  Nuclear  Physics 
Institute,  Nov.  1967. 
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Detection  of  Hydrogen  in  Semiconductors  by  Nuclear  Methods 

Thomas  Wichert 

Fakuhatfur  Physik,  Universiiax  Konstanz,  D-7750  Konstanz,  FRG 


The  impact  of  nuclear  methods,  like  the  ion  channeling  or  blocking  technique,  perturbed  ~n 
angular  correlation  (PAC)  and  Mossbauer  spectroscopy  (MS),  on  the  detection  of  hydrogen 
interacting  with  dopant  atoms  in  semiconductors  is  discussed.  In  particular,  the  recent  results 
for  H(D)-acceptor  interactions  in  Si  obtained  by  ion  channeling  and  PAC  are  be  presented 
and  their  interpretations  are  compared.  In  addition,  the  possibility  of  transferring  these 
experiments  to  studies  of  dopant  atom-hydrogen  interactions  in  GaAs  are  outlined. 

There  exist  several  ion  channeling  experiments  on  p-type  Si  doped  with  B  and  D  atoms  of 
equal  concentrations.  It  is  concluded  that  the  majority  of  D  atoms  is  situated  at  near-bond 
centre  sites  at  30  K  and  at  295  K;  however,  the  possibility  of  an  antibonding  site  of  D  at  295  K 
has  to  be  considered,  as  well.  The  interpretation  of  the  ion  channeling  results  is  affected  by 
two  experimental  facts:  The  high  D  (and  B)  concentration  (1019  cm"3),  required  for  this  type 
of  experiment,  and  the  introduction  of  the  D  atoms  from  a  plasma,  which  might  create 
additional  lattice  defects  being  capable  of  interacting  with  D  atoms,  as  well.  Recent  results  on 
the  lattice  location  of  radioactive  interstitial  impurities,  like  Li  atoms  in  GaAs,  using  the 
blocking  technique,  wfl]  be  presented  because  this  type  of  experiment  is  sensitive  to  much 
lower  impurity  concentrations  than  the  ion  channeling  technique. 

The  formation  of  acceptor-H  complexes,  involving  the  radioactive  probe  atom  niIn,  was 
identified  with  help  of  the  PAC  technique.  This  type  of  experiments  yields  information  on  the 
chemical  and  structural  identity  of  the  formed  complexes.  Using  the  H  specific  electric  field 
gradient,  the  formation  and  dissolution  of  the  In-H  complexes  was  pursued  under  a  variety  of 
different  conditions,  such  as  sample  temperature,  dopant  atom  concentration  and  position  of 
the  Fermi  level.  By  additionally  doping  the  samples  with  a  second  type  of  acceptor  atom  the 
stability  of  other  acceptor-H  complexes,  here  of  B-H  complexes,  was  determined.  In  order  to 
complement  the  investigations  on  H,  experiments  on  Li  atoms  and  an  interstitial  X  defect, 
produced  during  wafer  polishing,  were  performed.  The  X  defect  bears  close  resemblance  to 
the  H  defect  in  some  respects,  but  is  different  from  H  because  of  its  mobility  and  stability  at 
the  acceptor  atoms.  For  a  transfer  of  this  type  of  experiments  to  GaAs  radioactive  atoms  are 
needed  which  act  as  acceptor  or  donor  atoms  and  at  the  same  time  are  suited  for  either  PAC 
or  MS  spectroscopy.  In  case  of  PAC,  Cd  is  a  possible  candidate,  in  ease  of  MS,  Sn,  Sb,  and  Te 
can  serve  as  probe  atoms.  Thus,  experiments  on  H  in  GaAs,  being  similarly  detailed  as  in  Si, 
should  become  feasible  in  the  near  future. 
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]on  Beam  Studies  of  Hydrogen  in  Semiconductors 

Alan  D.  Marwick 

IBM  Research.  Division,  T.J.  Watson  Research  Laboratory, 

P.O.  Box  218,  Yorktown  Heights,  NY  10598 

Icn  beam  analysis  techniques,  including  nuclear  reaction  analysis,  channeling,  and 
Rutherford  backscattering  are  being  applied  to  several  aspects  of  the  behavior  of  hy¬ 
drogen  in  crystalline  semiconductors  and  simple  device  structures. 

Nuclear  reaction  profiling  of  hydrogen  in  MOS  capacitor  structures  with  .41  and  Au 
gates  was  used  to  show  that  most  of  the  hydrogen  in  such  structures  is  localized  in  a 
layer  at  the  gate-oxide  interface  [1].  This  hydrogen,  whose  amount  in  the  samples 
studied  varied  between  10“  and  3  x  10“  atoms/cm1,  is  probably  the  source  of  the 
hydrogenic  species  believed  to  diffuse  towards  the  oxide, 'Si  interface  during  hot  electron 
injection  or  after  exposure  to  ionizing  radiation.  The  ruesults  suggest  that  the  source  of 
the  H  in  the  samples  teas  air  exposure  of  the  oxide  before  metallization. 

The  lattice  location  of  hydrogen  (3/f)  and  boron  atoms  in  boron-hydrogen  complexes  in 
silicon  was  determined  using  ion  channeling  combined  with  nuclear  reaction  analysis. 
[2,  3],  The  boron  atom  was  found  to  be  slightly  displaced  off  a  lattice  site,  by 
0.28±0.03A,  while  the  majority  of  the  hydrogen  atoms  occupied  bond-center  positions. 
THe  samples  were  characterized  by  SIMS  and  Monte-Carlo  modelling  was  used  in  the 
analysis  of  the  channeling  data.  The  results  are  in  qualitative  agreement  with  theoretical 
models  of  the  structure  of  the  B-H  center. 

Present  work  includes  investigation  of  damage  induced  by  high  hydrogen  concentrations 
in  silicon,  channeling  studies  of  other  solute-hydrogen  complexes,  and  annealing  studies 
on  solute-hydrogen  complexes. 

[1]  A.D.  Marwick  and  D.R.  Young,  “Measurements  of  hydrogen  in  metal-oxide- 
semiconductor  structures  using  nuclear  reaction  profiling”,  J.  AppL  Phys.  63 
(1988),  2291. 

[2]  A  D.  Marwick,  G.S.  Oehrlein  and  N.M.  Johnson,  "Structure  of  the  boron- 
hydrogen  complex  in  crystalline  silicon",  Phys.  Rev,  B  36  (1987),  4539. 

[3]  A.D.  Marwick,  G.S.  Oehrlein,  J.H.  Barrett  and  N.M.  Johnson,  “The  structure  of 
the  boron-hydrogen  complex  in  silicon",  Mat.  Res.  Soc.  Symp.  Proc.  104  (1988), 
259. 
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C.  G.  Van  de  Walle  (Philips  Briarcliff  Manor):  "Theory  of  hydrogen  diffusion  and 
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caculations".  Invited  talk. 
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The  Structure  and  Properties  of  H  Complexes  in  GaAs.a-Si  and 

Diamond. 

P  Briddon,  R  Jones.  CMS  Lister 

Dept,  of  Physics,  University  of  Exeter.  Exeter.  EX44QL.LK 
Electronic  mail  jones(2'uk.ac.exeter 

We  hav  e  calculated  75.6.7/  the  structure  and  properties  of  H  in  clusters  of  atoms  simulating  the 
solid  state.  These  have  contained  50-90  atoms  and  have  surfaces  saturated  with  H.  We  use  a 
self-consistent  local  density  functional  method  with  s-  and  p-  Gaussian  orbitals  as  a  basis  for  the 
wave-functions  and  s-  Gaussian  orbitals  as  a  basis  for  a  fit  to  the  charge  density  / 1 .2.3.4/ .  Bond 
centred  orbitals  are  used  in  addition  to  atom  centred  ones. 

The  cluster  and  bases  are  sufficiently  large  that  bond  lengths  and  angles  are  accurately  found  to 
within  a  few  per-cent.  Dynamical  properties  are  given  to  within  10%.  The  cluster  size  is . 
however,  inadaquate  to  give  the  theoretical  band  gap,  but  we  feel  the  method  can  distinguish  deep 
from  shallow  levels. 

We  have  found  the  GaAs.Si  and  diamond  bond  lengths  to  be  within  25c,  1.59c.  15c  of  experiment 
respectively .  Calculating  the  energies  of  distorted  central  atoms  of  the  cluster  allowed  bond 
stretching  and  bending  force  constants  to  be  found  which  ,  in  the  case  of  Si  and  C,  yield  bulk 
phonon  dispersion  curves  within  10%  at  the  zone-edges.  For  GaAs ,  a  Coulomb  potential  must  be 
added  to  gi\’e  phonon  frequencies  within  20%  at  the  zone  edges. 

1.  H  in  GaAs,  GaAstBe,  GaAstP  /7/ 

56  or  57  atom  clusters  revealed  there  is  a  marginal  (.2eV)  tendency  for  neutral  H  to  lie  at  a 
T-intersdtia]  position  over  a  BC-(bond  centred)  site.  Furthur  calculations  on  larger  clusters  are 
necessary  to  confirm  the  stability  of  the  T-site.  BC  H  has  a  singly  occupied  level  high  up  in  the 
gap;  T-sited  H  has  a  level  in  the  middle  of  the  gap.  Thus  in  p-GaAs  ,  BC-  sited  H  is  likely  to  be 
preferred  over  T-sited  H  and  vice-vena  in  n-GaAs. 

Be/Ga  is  an  off-centred  defect  which  is  displaced  along  (-1,-1,-1)  and  bonded  to  3  As  atoms.  The 
As  atom  along  [1,1,1]  is  also  3-fold  coordinated  and  is  the  centre  for  the  acceptor  wave-function  . 
A  H  atom  prefers  to  sit  at  a  BC  site  between  the  As  and  Be  where  it’s  w-ave  function  strongly 
overlaps  that  of  the  acceptor. 

This  structure  and  the  LVM  of  H  (2083  cm'1)  are  in  good  agreement  with  experimental  values 
(2037  cm'1). 

Conversely  H  in  Si/Ga  prefers  to  lie  along  [-1,-1,-!]  where  the  Si  is  4-fold  coordinated.  The  Si-As 
bond  along  [1,1,1]  is  broken  and  the  low  lying  filled  state  has  appreciable  overlap  on  this  As  atom 
and  the  H.  The  calculated  H  stretch  mode  frequency  for  an  87  atom  cluster  is  1 805  cm- 1  .  The 
obsetx-edone  is  1717cm-l  . 

2.  H  in  a-Si  15/ 

A  34  atom  cluster  of  Si  and  H  atoms  containing  a  central  Si-H  bond  adjacent  to  a  strained  Si- Si 
bond  simulates  a-Si:H.  The  energy  of  this  cluster  was  evaluated  as  a  function  of  the  strain  in  the 
Si-Si  bond.  At  around  40%  strain  ,  the  Si-H  bond  spontaneously  breaks  and  the  H  atom  forms  a 
bond  centred  H  defect  in  the  middle  of  the  strained  Si-Si  bond  leaving  a  dangling  bond  behind.To 
move  the  H  atom  furthur  away  from  the  bond  -  centred  site  adjacent  to  the  Si  dangling  bond  costs 
-leV,  This  calculation  supports  a  model  of  the  Staebler-Wronski  effect  introduced  by  Stutzmann. 
Jackson  and  Tsai,  Light  of  energy  1 .2- 1 . 5e  V  creates  e-h  pairs  at  strained  bonds  w  hich  either 
promote  Si-H  dissociation  (the  first  step  above)  or  cause  a  furthur  Si-Si  weakening  allowing 
spontaneous  Si-H  dissociation  . 

3.  H  in  Diamond.  16/ 

A  57  atom  cluster  showed  BC  -sited  H  is  1 .9eV  more  stable  than  T-sited  H.  This  is  consistent 
with  muon  spin  rotation  experiments.  The  BC  H  has  a  level  close  to  Ec  and  thus  there  is  a  strong 
tendency  for  it  to  pair  with  acceptors  or  other  defects.  The  w  ave  function  is  built  up  of 
antibonding  p-orbitals  on  adjacent  C  atoms  and  has  lobes  which  point  along  [-1-1-1]  and 
consequently  there  is  strong  overlap  w  it h  H  placed  at  an  anti-bonding  sue.  This  structure  of  H2  is 
more  stable  by  3.32  eV  over  a  Hj  molecule  placed  near  a  T  site.  Accordingly  H2  should  should 

have  a  deep  trapped  state  in  diamond. 
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Theory  of  hydrogen  diffusion  and  reactions 
in  crystalline  semiconductors 

Chris  G.  Van  de  Walle 
Philips  Laboratories 
Briarcliff  Manor,  New  York 

State-of-the-art  theoretical  techniques  can  provide  a  first-principles  description 
of  the  behavior  of  hydrogen  in  crystalline  semiconductors.  Using  density-functional 
theory  and  ab  initio  pseudopotentials,  a  comprehensive  picture  can  be  obtained  of 
stable  configurations,  migration  paths,  charge  states,  etc.,  for  H  in  pure  and  in  doped 
semiconductors.  Hydrogen  interacts  strongly  with  the  crystalline  network,  and  re¬ 
laxation  of  the  lattice  around  the  impurity  must  be  included  to  obtain  the  correct 
minimum-energy  positions.  Total-energy  surfaces,  in  which  the  position  of  the  H 
atom  in  the  crystal  is  the  coordinate,  provide  immediate  insight  into  the  relevant 
properties.  Silicon,  as  the  prototypal  semiconductor,  has  been  the  first  material  for 
which  we  have  performed  the  calculations1;  recent  work  on  H  in  Ge2  indicates  similar 
behavior. 

The  study  of  H  as  an  impurity  in  ihe  semiconductor  provides  information  about 
its  charge  state,  and  allows  us  to  predict  the  reaction  mechanisms  through  which  it 
will  interact  with  other  impurities.  For  instance,  H  in  p-type  Si  is  most  stable  in 
the  positive  charge  state;  its  electron  causes  compensation  of  free  holes.  Such  an 
H+  impurity  diffuses  very  fast  through  the  crystal,  and  will  readily  be  attracted  by 
negatively  charged  acceptors.  Formation  of  neutral  hydrogen-acceptor  pairs  is  a  result 
of  this  process.  Our  calculations  of  boron-hvdrogen  complexes  in  Si  show  that  the 
global  energy  minimum  occurs  for  H  at  the  center  of  a  Si-B  bond.  H  can  move  on 
a  spherical  shell  around  B  with  an  energy  barrier  of  only  0.2  eV.  In  n-type  Si,  we 
predict  H  to  be  negatively  charged.  The  deep  doDor  and  acceptor  levels  associated 
with  the  amphoteric  H  impurity  are  close  in  energy;  the  possibility  that  H  in  Si  is  a 
negative-U  impurity  can  not  be  excluded. 

In  Si,  isolated  hydrogen  atoms  in  the  positive  and  neutral  charge  state  prefer  the 
bond-center  site  as  the  stable  location;  in  the  negative  charge  state,  the  tetrahedral 
interstitial  site  is  favored.  In  the  absence  of  other  defects,  H2  molecules  are  the  most 
stable  state  for  H  in  the  Si  crystal.  The  binding  energy  (compared  to  two  >'  W 
neutral  H  atoms  inside  the  crystal)  is  2  eV  per  molecule.  If  one  of  the  H  atoms  that 
forms  the  molecule  is  positively  charged  (as  suggested  by  Johnson  and  Herring3)  the 
binding  energy  becomes  smaller,  and  is  Fermi- level  dependent. 


1  C.  G.  Van  de  Halle,  Y.  Bar- Yam.  and  S.  T.  Pantelides,  Pbvs.  Rev.  Lett.  60.  2761 
(1988). 

2  P.  J.  H.  Denteneer.  C.  G.  Van  de  Halle,  and  S.  T.  Pantelides.  to  be  published. 

3  N.  M.  Johnson  and  C.  Herring,  Phys.  Rev.  B  38.  1581  (19S8). 
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Theory  of  Interstitial  Hydrogen  and  Hydrogen-Containing 
Aggregates  in  Silicon:  Recent  Semiempirical  Calculations* 

Gary  G.  DeLeo 

Department  of  Physics  and  Sherman  Fairchild  Center 
Physics  Building  #  16 
Lehigh  University 

Bethlehem,  Pennsylvania  18015  USA 


The  results  of  recent  seaiempirical  computations  on 
hydrogen-related  defects  in  silicon  are  summarized.  Total 
energies  and  single-particle  electronic  structures  are  produced 
by  the  KNDO  (modified  neglect  of  diatomic  overlap)  method  using  a 
finite-cluster  framework.  Equilibrium  geometries,  migration  and 
reorientation  barriers,  and  vibrational  frequencies  are  inferred 
from  the  atom-position  dependence  of  the  total  energies. 
Specifically,  interstitial  hydrogen,  hydrogen- acceptor  (B,  Al) 
pairs,  and  a  hydrogen-donor  (P)  pair  are  treated. 

The  impurities  and  environment  are  simulated  in  these 
treatments  by  finite  fragments  of  the  solid,  or  "clusters."  The 
primary  cluster  used  in  these  studies  is  XSi  8H18 '  this  is  bond 
centered.  Also,  we  use  the  clusters  XSi10H1$  and  XSi4H12  to 
simulate  interstitial  and  substitutional  environments, 
respectively.  The  hydrogens  serve  to  terminate  the  cluster. 

In  its  simplest  application,  the  host  atoms  are  kept  at 
their  rigid-lattice  sites  and  the  total  energy  is  monitored  as  a 
function  of  the  impurity  positions.  A  somewhat  more  elaborate 
procedure  involves  the  simultaneous  displacments  of  the  impurity 
atom(s)  and  the  neighboring  silicon(s) .  When  the  outer  silicons 
are  relaxed,  they  are  connected  to  springs  with  parametrized 
spring  constants  for  parallel  (to  bond)  and  perpendicular  (to 
bond)  displacements;  we  refer  to  this  as  the  hybrid  approach. 

The  H-terminators  follow  the  outer  silicons  with  the 
corresponding  bond  angles  and  lengths  preserved  (taken  here  to  be 
109.47°  and  1.48  A,  respectively).  We  note  that  allowing 
silicons  to  fully  relax  for  each  hydrogen  position  may  not  be 
appropriate  when  considering  the  dynamics  of  a  light  element  such 
as  hydrogen. 

The  most  decisive  statement  which  can  be  made  on  the  basis 
of  our  computational  results  for  isolated  interstitial  hydrogen 
is  that  a  bond-centered  (BC)  configuration  is  the  stable  geometry 
for  the  neutral  and  positive  charge  states. [1]  This  supports 
other  computational  results  [2]  and  recent  >^SR  [3]  and  EPR 
studies  of  anomalous  muonium  and  hydrogen  in  silicon, 
respectively.  We  find  neutral  hydrogen  to  be  asymmetrically 
located  along  the  bond,  with  a  barrier  for  reorientation  of  about 
0.2  ev.  This  barrier  is  sufficiently  small,  however,  that 
motional  averaging  would  cloud  an  inequivalence  of  the 
neighboring  two  silicons.  Such  an  asymmetric  configuration  might 
be  distinguished  from  the  corresponding  symmetric  one  in  an 
infrared  experiment  by  the  anharmonic  H/D  isotope  shift  of  1.7 
which  we  predict. 
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Since  the  neighboring  silicons  relax  by  about  0.4  A  each  to 
accommodate  the  hydrogen,  migration  between  neighboring  BC  sites 
would  require  substantial  silicon  participation.  Therefore,  we 
compute  a  migration  barrier  by  determining  the  energy  required  to 
displace  the  neighboring  silicons  so  that  the  hydrogen  may 
spontaneously  move  between  adjacent  BC  sites.  We  find  the 
migration  barrier  to  be  ^1.0  eV. 

There  also  appears  to  be  an  athermal  motion  which  is 
associated  with  the  migration  of  normal  muonium.  Here  it  is  most 
appropriate  to  determine  the  migration  barrier  in  the  absence  of 
silicon  displacement;  the  hydrogen  would  tunnel  before  the 
silicons  could  fully  respond. [4]  We  find  a  hydrogen-motion 
barrier  for  migration  between  adjacent  M-sites  of  ~0.1  -  0.2  eV. 
Another  way  of  expressing  such  athermal  motion  may  involve  a 
Bloch  representation  of  the  H  or  muonium  atoms;  this  is  currently 
being  explored. 

In  some  recent  treatments  of  the  hydrogen-acceptor  pairs,  we 
have  taken  up  the  problem  of  the  low  frequency  (^70  cm-1) 
sideband,  which  was  revealed  experimentally  by  Stavola,  et  al.[5] 
We  have  considered  the  hydrogen  off-axis  motion  (wagging,  or  E 
mode)  as  a  candidate  for  this  sideband.  We  compute  these 
vibrational  frequencies  to  be  in  the  700-800  cm-1  range,  much  too 
high  to  explain  the  observed  sideband.  Furthermore,  they  appear 
to  be  relatively  insensitive  to  the  locations  of  the  near¬ 
neighbor  silicons  and  whether  or  not  the  hydrogen  is  forced  off 
the  axis  by  reducing  the  Si-acceptor  distance.  Hence,  it  is  not 
clear  that  these  wagging  modes  can  explain  the  low-frequency 
sideband.  Also,  more  complicated  motions  involving  the  hydrogen 
are  not  likely  to  exhibit  an  isotope  shift  so  close  to  1.41. 

This  remains  a  problem  of  interest  and  controversy. 

In  some  very  preliminary  work,  we  have  explored  the  energy 
surface  appropriate  to  hydrogen  motion  in  a  hydrogen-phosphorous 
pair.  Current  models  locate  the  hydrogen  along  the  <lll>  axis  in 
an  antibonding  position  to  the  silicon  neighboring  the 
phosphorous. [6]  We  have  seen  indications  that  the  hydrogen  may 
be  stable  as  well  at  a  bond-centered  site.  These  studies  are 
continuing  with  larger  clusters. 

*  Work  done  in  collaboration  with  W.  Beall  Fowler. 
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A  tight  binding  total  energy  calculation  indicate  a  high 
chemical  activity  cf  hydrogen  ; n  a  Si  lattice.  Atomic  hydrogen 
is  not  stable  m  an  unbound  interstitial  form  m  a  Si  lattice 
and  reacts  with  the  Si— Si  bond,  as  the  result  of  such  reaction 
an  H-center  is  created  whitch  consist  a  o  type  Si(p)-H(s)  bond 
and  one  unsaturated  p-like  dangling  hybrid.  The  hydrogen  atom 
may  be  bound  outside  or  inside  the  bond  region  depending  on 
the  occupancy  of  the  adjacent  dangling  hybrid.  The  calculated 
activation  energy  of  the  diffusion  cf  H-centers  is  1.  £3  eV. 
Interactions  between  H-centers  and  D-centers  lead  to 
annihilation  of  dangling  hybrids  m  the  lattice,  in  spite  of 
the  fact  that  the  number  of  Si-H  bonds  remains  unchanged.  On 
the  o+her  hand,  collisions  between  diffusing  H-centers  may 
lead  to  two  types  of  effects:  1)  formation  of  an  Hg  molecule 
in  the  Si  lattice,  or  £)  formation  of  a  stable  defect  contain 
two  Si ( p ) -H ( s )  bonds.  The  first  of  the  predicted  effects  gives 
fresh  microscopic  insight  .  into  the  mechanism  of  lew 
temperature  evolution  of  hydrogen  from  a-Si:H.  The  second 
effect  is  responsible  for  the  damage  induced  by  hydrogen  m 
crystalline  Si. 
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